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FOREWORD 


Cornell  Aeronautical  Laboratory,  Inc.  (CAL),  Buffalo,  New  York, 
performed  che  research  reported  herein  for  Air  Programs,  Naval  Applications 
Group,  Office  of  Naval  Research,  Washington,  D.  C.  20360,  under  ONR 
Contract  Nonr  3691(00).  Dr.  P.  Crimi  of  CAL  conducted  the  study,  with 
assistance  from  Mr.  W.  G.  Brady.  Mr.  Harvey  Selib  was  responsible 
for  the  preparation  of  the  digital  machine  program.  The  program  was  under 
the  supervision  of  Mr.  F.  DuWaldt  who  prepared  this  summary  report. 

The  work  was  performed  from  February  1965  through  May  1966 
and  is  a  continuation  of  the  work  described  in  CAL  Report  BB-1665-S-2, 
"Representation  of  Propeller  Wakes  by  Systems  of  Finite  Core  Vortices", 
by  W.  G.  Brady  and  P.  Crimi,  issued  under  ONR  Contract  Nonr  3691(00), 
February  1965  (DDC  No.  612007).  As  an  outgrowth  of  this  research,  a 
companion  study  on  rotor  flow  lields  (in  the  absence  of  a  ground  effect)  was 
performed  for  the  U.  S.  Army  Ballistic  Research  Laboratories  (Contract 
No.  DA  30-069-AMC-645(R))  by  Dr.  P.  Crimi  and  Mr.  A.  Sowyrda  and 
published  in  two  volumes: 

1.  Crimi,  P.  Theoretical  Prediction  of  the  Flow  in  the 
Wake  of  a  Helicopter  Rotor  Part  1  -  Development  of 
Theory  and  Results  of  Computations  CAL  Report 
No.  BB-1994-S-1  September  1965 

2.  Crimi,  P.  Theoretical  Prediction  of  the  Flow  in  the 
Wake  of  a  Helicopt^  Rotor  Part  2  -  Formulation  and 
Application  of  the  Rotor -Wake  Flow  Computer  Program 
CAL  Report  No.  BB-1994-S-2  September  1965 


m 


BB-1665-S-3 


The  analyses  reported  herein  (in  which  the  effect  of  ground  proximity  is 
included)  are  similar  to  those  used  in  the  companion  study  since  both  stem 
from  the  model  developed  by  W.  G.  Brady,  P.  Crimi,  and  A,  Sowyrda 
under  ONR  Contract  Nonr  3691(00). 
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ABSTRACT 


This  is  a  summary  report  that  presents  the  development  of  a  wake 
model  for  a  lifting  propeller  (rotor)  in  ground  effect  and  the  computational 
procedure  used  to  determine  the  spacial  distribution  of  wake  vorticity  and 
the  induced  velocity  field  accompanying  that  vorticity  distribution. 

Sample  calculations  for  a  two-bladed  rotor  were  carried  out  on  an 

IBM  7044  computer.  Locations  of  wake  vortical  elements  and  the  associated 

induced  velocities  at  selected  field  points  are  presented  for  advance  ratios 

H 

of  0,  0.  02,  0.  05,  and  0.  10  for  an--  (ratio  of  rotor  height  above  ground  tc. 

H 

rotor  radius)  of  1.0.  Also,  a  few  results  for  a  hovering  case  with  —  =  0.  5 

.K 

are  presented.  A  calculated  root-mean-square  velocity  map  is  compared 
with  measured  hovering  data  (time  averaged)  and  good  agreement  is  obtained 
’n  the  outer  half  of  tho  slipstream.  The  implications  of  computed  results 
with  respect  to  the  flov  field  at  a  horizontal  stabilizer  or  tail  rotor  are  noted 
as  examples  of  the  use  ol  the  model. 

The  wake  vortex  model  used  for  these  calculations  is  a  distorted 
continuous  helix  emanating  from  each  blade  tip.  For  purposes  of  numerical 
calculation  of  mutual  vortex  interference  effects,  the  continuous  vor^x  is 
approximated  by  short  straight-line  filaments.  For  purposes  of  numerical 
calculation  of  self-induced  effects  (that  is,  the  effect  of  an  element  on  itself), 
the  local  element  curvature  and  the  core  diameter  are  estimated.  The 
Fortran  statements  of  the  computer  program  are  presented  in  the  A.ppendix. 
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SYMBOLS 


Of 

ratio  of  vortex  core  radius  to  rotor  radius,  — 

R 

H/R 

ratio  of  height  of  rotor  hub  above  ground  to  rotor  radius 

N 'b 

number  of  rotor  blades 

p 

air  density- 

p 

designation  for  ft  wake  reference  point 

R 

rotor  radius 

& 

radius  of  curvature  of  a  vortex  tube 

r 

position  vector  of  a  point  in  the  flow 

At 

incremental  time 

self -induced  velocity  of  a  vortex  element 

components  of  fluid  velocity  in  the  directions  of  x>,l/  and 
respectively,  nondimensionalized  by  X  (1R 

^  'i/  »  K/t-  '  Vlt) 

time  averages  of  1/ ;  l/,  and  respectively 

It 

free-stream  velocity 

\L 

local  fluid  velocity 

Vr 

time  average  of  the  radial  velocity 

if 

induced  downwash 

~R 

time  average  of  if 

% 

momentum  value  for  the  time  average  of  if 

W 

aircraft  weight 

(x.y,}) 

rectilinear  coordinates  with  origin  in  the  tip-path  plane, 

nondimensionalized  by  R 

Xll 


BB-1 665-S-3 


a,T 

r 

A 

M' 

V 

a 


inclination  of  the  free  stream  to  the  tip-path  plane 
circulation 

blade  loading  parameter;  A  =  ‘f'Wf (tTZN g yOCi2 R*) 
advance  ratio;  -  Vf/(RtR) 
rotor  azimuth  angle 
rotor  angular  speed 

magnitude  of  distributed  trailing  vorticity  in  a  wing  wake 
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INTRODUCTION 


Position  and  structure  of  the  vortical  wake  of  a  propeller  (rotor) 
are  important  considerations  in  many  practical  problems.  Among  these 
are  (1)  calculation  of  rotor  aerodynamics  for  the  purpose  of  estimating 
performance,  blade  dynamic  response,  and  stability  and  control  derivatives; 
(2)  prediction  of  the  flow  field  in  which  auxiliary  devices  operate;  (3)  esti¬ 
mation  of  rotor -induced  loads  on  aircraft  structures  other  than  the  rotor; 

(4)  estimation  of  the  potential  for  entrainment  and  transport  of  ground  debris 
or  water  droplets;  and  (5)  correction  of  wind-tunnel  data  for  wall  effects. 

When  an  actuator  is  in  an  axial-flight  condition  and  the  free-stream  velocity 
is  large  compared  to  ihe  wake-induced  velocities,  approximate  techniques 
can  be  used  in  analyzing  the  effect  of  the  wake  structure  relative  to  the 
various  problem  areas  indicated  above.  However,  when  an  actuator  is  in 
slow  speed  flight  and/or  is  highly  yawed,  the  wake  structure  becomes  an 
important  influence.  Of  particular  importance  in  this  latter  category  are 
the  flow  fields  associated  with  VTOL  aircraft  in  the  near-hovering  condition 
when  the  actuator  axis  is  approximately  perpendicular  to  a  nearby  ground 
plane.  It  is  to  this  problem  that  the  current  research  was  directed. 

Mathematical  description  of  the  flow  field  for  a  rotor  is  made  diffi¬ 
cult  by  the  interdependence  of  the  solution  and  the  wake  position  --  that  is, 
the  boundary  condition  is  known  but  the  location  of  the  boundary  is  not  known. 
This  situation  is  characteristic  of  problems  that  can  be  generally  classified  as 
"free -boundary  problems".  The  method  of  solution  that  can  sometimes  be 
used  on  this  class  of  problems  is  to  guess  the  position  of  the  boundary, 
compute  the  solution,  and  determine  if  the  computed  solution  is  consistent 
with  the  assumed  boundary  location.  This  approach  is  based  on  the  use  of 
space  coordinates  as  independent  variables  and  the  computation  is  cast  as 
some  sort  of  feedback  scheme  in  which  the  differences  between  the  assumed 
boundary  location  and  the  resulting  solution  are  used  to  determine  a  new 
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estimated  boundary  position.  The  technique  involved  in  the  current  studies 
is  slightly  different.  Advantage  is  taken  of  the  azimuthal  periodicity  of  the 
rotor  position  and  the  independent  variable  is  chosen  m  be  time  (or  azimuthal 
angle  of  a  reference  blade).  The  point  of  view,  then,  is  that  "new"  elements 
of  wake  vorticity  are  generated  as  the  blading  rotates  and  translates.  These 
elements  are  convected  with  velocities  determined  by  their  self-induced 
velocities  as  well  as  velocities  induced  by  the  existing  wake  structure,  the 
bound  vorticity  of  the  blade,  and  the  image  vortices  introduced  to  satisfy  the 
ground-plane  boundary  condition.  It  is  reasonable  to  assume  that,  eventually, 
a  stabilized  periodic  wake  array  would  be  obtained  since  the  condition  of  a 
fixed  periodic  blade  loading  was  imposed.  This  is,  however,  not  always  the 
case  --  this  situation  will  be  pointed  out  in  the  discussion  of  the  sample  cal¬ 
culations.  It  is  believed  that  the  lack  of  perfect  periodicity  in  the  solutions 
is  a  result  of  the  truncation  of  the  wake  extent;  but  the  possibility  of  a  compu¬ 
tational  --or  even  physical  --  instability  for  some  cases  cannot  be  completely 
discounted. 

Initial  effort  in  this  research  task  was  devoted  to  the  special  case 
of  actuators  hovering  near  a  ground  plane  (Reference  1).  The  wake  structure 
was  assumed  to  be  representable  by  discrete  concentric  planar  vortex  rings 
of  finite  core  diameter.  Image  rings  were  used  to  enforce  the  ground-plane 
boundary  condition.  This  model  subsequently  has  been  used  in  the  estimation 
of  slipstream  envelopes  for  propellers  at  zero  advance  ratio  (Reference  2). 

Vortex  rings  were  useful  devices  in  analyzing  hovering  cases 
because  the  requirement  for  axial  symmetry  gave  assurance  that  ihe  rings 
would  remain  planar.  In  the  presence  of  a  free-stream  velocity  having  a 
component  parallel  to  the  actuator  plane,  however,  the  axial  symmetry 
disappears  and  the  advantage  of  the  ring  representation  is  lost.  Therefore, 
the  continuous  distorted  helix  (recommended  in  Reference  1)  was  selected  as 
the  model  for  the  current  effort.  Again,  the  ground-plane  boundary  condition 
was  enforced  by  the  influence  of  an  image  wake. 


2 


BB-1 665-S-3 


A  recently  completed  study  (References  3  and  4)  has  shown  that 
the  distorted  helix  model  yields  time-averaged  velocities  that  agree  well 
with  measurements  made  in  and  about  a  helicopter  rotor  wake  (in  the  absence 
of  an  appreciable  ground  effect).  No  check  could  be  made  of  the  ability  of 
the  model  to  represent  the  time  variations  of  the  velocities  because  measure¬ 
ments  are  not  available. 

The  digital  machine  computing  program  is  based  on  the  assumption 
that  the  rotor  is  in  steady  level  flight  (or  hovering)  with  specified  tip-path- 
plane  angle.  Shaft  rotational  speed,  rotor  force,  initial  vortex  core  diam¬ 
eter,  number  of  blades  and  the  ratio  of  rotor  radius  to  height  above  the 
ground  are  also  inputs.  Information  relative  to  the  computing  program  is 
given  in  the  Appendix. 
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GENESIS  OF  THE  IDEALIZED  MODEL 


A  detailed  analysis  of  a  propeller  or  rotor  flow  field  in  general 
form  would  be  prohibitively  expensive  in  terms  of  man-power  required 
and  computing  costs.  It  is  necessary,  therefore,  to  idealize  the  physical 
situation  in  order  to  reduce  the  mathematical  complexity.  Ultimate  justi¬ 
fication  for  the  validity  of  the  simplifications  introduced  rests  upon  the 
comparison  of  the  essential  quantities  with  experimental  observables  or 
comparison  with  a  more  representative  analytical  model. 

The  fluid  velocities  depend  primarily  upon  the  effects  of  the  rotor 
blades,  their  wake,  and  the  ground  boundary  proximity.  The  effects  are,  of 
course,  interrelated  but,  for  purposes  of  discussion,  will  be  considered 
separately. 

Rotor  or  propeller  blades  may  be  regarded  as  high -aspect-ratio 
wings  in  a  harmonically  varying  stream.  The  actuator  lift  must  be  approxi¬ 
mately  equal  to  the  weight  of  the  aircraft.  Therefore,  the  circulation  about 
a  given  blade  must  be  a  harmonic  function  of  azimuth.  The  radial  variation 
of  circulation  depends  on  blade  twist,  planform,  and  the  influence  of  the 
wake.  Nonetheless,  results  of  theoretical  computations  of  harmonic  air 
loads  and  measurements  (Reference  5)  indicate  that  the  radial  variation 
of  circulation  is  approximately  elliptical. 

The  wake  of  each  blade  is  a  thin  sheet  of  vortical  fluid.  When  the 
rotor  is  translating,  the  vorticity  in  the  sheet  has  two  components  (shed 
vortices  initially  parallel  to  the  spanwise  direction  and  trailing  vortices 
initially  parallel  to  the  chordwise  direction).  Flow  visualization  studies 
indicate  that  the  sheet  rolls  up  within  a  few  chord  lengths  of  the  blade  to 
form  a  single  concentrated  vortex  trailing  from  near  the  tip  of  the  blade 
(see,  for  example,  Reference  6).  Theoreticall/,  a  similar  vortex  should 
arise  from  the  vicinity  of  the  blade  root.  However,  published  smoke  pictures 
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as  well  as  an  analytical  treatment  of  the  wake  of  a  hovering  rotor  lead  to 
the  tentative  conclusion  {Reference  1)  that  the  root  vortex  is  rapidly  con- 
vected  away  from  the  rotor  and/or  dissipated. 

The  effect  of  the  ground  can  be  accounted  for  by  the  introduction  of 
an  image  vortex  for  every  element  of  the  vorticity  in  the  fluid.  In  this  way, 
the  ground  boundary  condition  (the  vanishing  of  the  normal  velocity  com¬ 
ponent)  is  satisfied. 

The  above  observations  lead  to  the  idealized  model  used  for  compu¬ 
tational  purposes.  The  fluid  is  considered  to  be  inviscid  and  incompressible. 

First,  the  rotor  blades  are  replaced  by  radial  line  vortices  whose 
strengths  vary  once  per  revolution  (  T ~ /- 2 stn  ^  ),  Blade-induced  effects 
will  be  reasonably  reproduced  except  in  the  immediate  vicinity  of  the  blades. 
The  neglect  of  higher  harmonic  variations  of  blade  circulation  is  not  expected 
to  be  important  in  comparison  with  the  wake  effects  for  field  points  not  too 
close  to  the  blading.  As  a  consequence  of  this  representation,  the  blade 
loading  can  be  regarded  as  a  known  quantity. 

Second,  the  wake  is  idealized  as  individual  free  vortices,  one  trailing 
from  the  tip  of  each  blade.  The  circulation  of  each  tip  vortex  varies  in  accord¬ 
ance  with  the  (known)  azimuthal  variation  of  the  blade  circulation.  The  trailing 
vortices  are  assumed  to  have  circular  cross  sections  for  purposes  of  the 
computation  of  the  self-induced  velocity.  The  effects  of  shed  vortices  and 
trailing  root  vortices  (required  by  Kelvin's  theorem)  have  been  neglected  as 
noted  above. 

Corresponding  to  each  blade  bound  vortex  and  each  wake  vortex 
element,  there  is  an  image  vortex  below  the  ground  plane. 

Formulation  of  the  digital  model  for  the  flow  corresponding  to  the 
simplified  model  can  be  accomplished  in  a  straightforward  manner.  A 
coordinate  system  fixed  in  the  tip  path  plane  is  introduced  (see  Figure  1). 
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The  vector  V_p  in  Figure  1  is  the  free -stream  velocity  (the  negative  of  the 
translational  velocity  of  the  aircraft),  inclined  at  an  angle  ocT  to  the  x  -axis 
and  parallel  to  the  plane.  The  angle  denotes  the  azimuth  position  of 
the  reference  blade,  and  jr  is  the  vector  positioning  a  given  point  with  re¬ 
spect  to  the  origin. 

The  fluid  velocity  at  a  given  point  located  by  the  vector  r  may 
be  expressed  in  the  form 


r  (r)rf  xdr 


(1) 


where  r1  -rp-jr  and  r  is  the  circulation  about  the  vortex  element  at_r  . 
The  line  integral  is  to  be  taken  over  all  vortices  in  the  flow  ana  the  image 
system  --  which  paths  are  collectively  denoted  by  Cv  .  It  is  necessary  to 
modify  this  expression  when  =  £  and  this  point  is  discussed  later.  The 
circulation,  P ,  about  the  blade  vortices  is  estimated  directly  in  terms  of 
flight  parameters.  The  circulation  about  a  wake  vortex  at  a  given  point  is 
simply  that  value  assigned  to  the  blade  vortex  when  it  was  generating  that 
segment  of  wake. 

The  only  information  lacking  for  the  complete  specification  of  the 
flow  at  a  given  instant,  then,  is  the  location  of  the  wake  vortices  at  that 
instant.  The  position  of  a  given  point  on  a  wake  vortex  located  by  the 
vector  r  is  the  time  integral  of  the  velocity  experienced  by  that  fluid 
particle 


(2) 
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Thus,  even  after  being  simplified,  the  flow  can  only  be  obtained  as  the 
solution  of  the  nonlinear  integral  equation  formed  by  the  substitution  of 
Equation  (1)  into  Equation  (2).  A  direct  analytical  solution  is  not  feasible 
but  the  problem  is  amenable  to  solution  by  numerical  methods  using  a  high¬ 
speed  digital  computer.  The  manner  in  which  the  formulations  of  Equations 
(1)  and  (2)  were  implemented  for  digital  computation  is  discussed  in  the 
following  section. 
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DESCRIPTION  OF  THE  COMPUTATION  MODEL 


The  idealized  continuous  model  deduced  in  the  previous  section 
is  replaced  by  a  sectioned  model  having  locally  constant  characteristics. 
Consistent  with  this  replacement  is  the  representation  of  integrals  by  sums. 
This  same  technique  was  used  in  Reference  3  and  the  following  discussion 
follows  closely  that  given  in  Reference  3. 

Wake  Positioning  and  Displacement 

For  purposes  of  numerical  analysis,  a  series  of  points  is  marked 
off  on  each  wake  vortex,  as  shown  schematically  in  Figure  2. 


Figure  2  WAKE  REFERENCE  POINT  IDENTIFICATION 
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Points  are  placed  sufficiently  close  together  so  that,  for  purposes  of 
computation  of  wake-induced  velocities,  the  segment  (L)  of  wake  between 
points  Pi  and  may  be  considered  to  be  a  rectilinear  vortex  having  con¬ 
stant  circulation  along  its  length.  The  wake  spacial  configuration  at  any 
instant  is  then  defined  by  the  locations  of  the  end  points. 

Computations  are  initiated  by  specifying  that  each  of  the  wake 
vortices  lies  on  a  prescribed  curve.  Generally,  the  curve  is  chosen  to 
be  a  skewed  helix  with  pitch  and  skew  selected  from  simple  momentum 
considerations.  The  numerical  analogues  to  Equations  (1)  and  (2)  are 
then  performed  by  first  summing  the  velocity  contributions  of  all  vortex 
elements  in  the  flow  at  each  reference  point  (Equation  (1)),  and  then  con¬ 
vening  each  point  at  the  computed  velocity  for  a  time  interval  At  (Equation 
(2)).  The  time  interval  At  is  chosen  to  correspond  to  a  small  finite  change 
in  the  azimuth  position,  AR ,  of  the  blades 


At  « 


At> 

SI 


(3) 


where  fl  is  rotor  angular  speed.  Once  the  new  coordinates  of  each  reference 
point  are  computed,  the  azimuth  of  the  blade  vortices  is  increased  by  A/'  and 
the  velocity  computation  is  performed  again.  As  the  blade  vortices  are  re¬ 
positioned,  a  wake  vortex  element  is  added  to  the  flow  at  the  tip  of  each  blade 
vortex,  the  added  vortex  having  a  length  of  approximately  RAW,  where  R  is 
rotor  radius.  A  corresponding  element  is  dropped  from  the  computations  at 
the  downstream  end  of  each  wake  vortex,  so  the  program  is  not  encumbered 
by  a  wake  of  ever-increasing  size.  The  computations  are  continued  in  this 
manner  for  a  specified  time  and  the  results  are  inspected.  If  a  neaiiy  peri¬ 
odic  solution  is  established  in  the  space  volume  of  interest,  the  calculation 
is  terminated. 
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The  total  number  of  wake  vortices  taken  into  account  and  the 
magnitude  of  AP  determine  the  accuracy  of  the  flow  representation  at  a 
given  point.  It  is  judged  that,  for  a  two-bladed  rotor,  a  value  for  of 
thirty  degrees  is  sufficiently  small  to  furnish  an  acceptable  estimate  of 
the  time  variations  of  the  flow  consistent  with  the  other  approximations 
introduced.  The  number  of  wake  elements  considered  depends  on  the 
region  of  interest  and  the  forward  speed.  If  the  flow  "near"  the  rotor 
plane  is  desired,  then  the  wake  generated  by  only  a  few  revolutions  of 
the  rotor  is  needed.  It  has  been  found,  for  example,  that  with  an  advance 
ratio,  pc,  of  about  0,  25  just  two  revolutions  of  wake  sufficiently  define  the 
flow  near  the  rotor  plane  in  the  absence  of  a  ground  plane  (Reference  3). 

At  lower  advance  ratios,  the  free  stream  does  not  clear  the  wake  from 
under  the  rotor  as  rapidly  so  more  wake  must  be  retained. 

Computation  of  Velocity  Induced  by  Wake  and  Blades 

The  velocity  induced  at  an  arbitrary  point  P  by  the  vortices  repre¬ 
senting  the  wake  and  blades  is  simply  the  sum  of  the  effects  of  an  array  of 
rectilinear  vortex  segments.  If  V  denotes  the  velocity  induced  at  P  by 
the  element  between  points  and  P^i_j ,  it  is  found  from  Equation  (1)  that 


( cos  6;  -  cos  Oi+t)  (4) 

where  /}  is  the  strength  of  the  element  and  $,■  ,  and  h  are  defined  in 

Figure  3.  The  velocity  is  directed  normal  to  the  plane  containing  Pt ,  Pu+1 
and  P  . 


As  the  field  point  P  is  made  to  approach  any  point  on  the  line 
joining  and  in  Figure  3,  the  induced  velocity  increases  without 
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WAKE  VORTEX  ELEMENT  [l  ) 


Figure  3  GEOMETRIC  RELATIONSHIPS  DEFINING  THE  FLOW 
INDUCED  BY  A  RECTILINEAR  VORTEX  ELEMENT 


limit  because  h  tends  to  zero  (Equation  4);  the  velocity  becomes  indeter¬ 
minate  for  h  =  0.  The  possibility  that  this  situation  will  arise  at  an  arbi¬ 
trarily  chosen  field  point  is  remote.  Conversely,  precisely  this  situation 
arises  in  the  evaluation  of  the  transport  velocities  of  the  ends  of  each  of  the 
wake  vorticity  elements  --  that  is,  at  points  P-  and  Pi+i  in  Figure  3.  Now 
the  singular  behavior  can  be  considered  to  have  been  introduced  as  a  result 
of  a  limiting  process  in  which  a  vortex  element  having  a  finite  core  diameter 
retained  its  strength  (circulation)  as  the  core  diameter  was  allowed  to  ap¬ 
proach  zero.  Evidently,  this  limiting  process  is  not  an  admissible  operation 
when  the  self-induced  ■'■elocity  of  a  vortex  is  sought.  Therefore,  for  purposes 
of  computing  wake -induced  flow  at  a  given  wake  reference  point  --  i.  e.  ,  a 
vortex  joining  point  --  the  effect  of  the  two  wake  elements  adjacent  to  that 
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point  must  be  analyzed  by  taking  core  size  into  account  (if  a  point  lies  outside 
the  core  substance  of  a  vortex,  the  formula  for  induced  velocity  is  the  same 
as  for  a  vortex  with  infitesimal  core  size;  i.  e.  ,  Equation  (4)).  The  self- 
induced  velocities  also  depend  directly  on  the  curvature  of  (the  centerline 
of)  the  vortex  tube  at  the  point  in  question.  An  approximation  to  the  actual 
curvature  is  obtained  by  assuming  that  the  '.wo  elements  whose  effects  at 
point  P-  are  to  be  computed  lie  on  the  circular  arc  which  passes  through 
points  Pc  f  ,  P-  and  P-L  +  } ,  as  shown  in  Figure  4.  Then  it  is  assumed  that 
the  two  vortex  segments  have  circular  cores  of  rotational  fluid  of  radii  a-_t 
and  a.-,  respectively.  Further,  the  approximation  is  made  that  the  vorticity 
within  the  core  is  tangent  to  the  arc  of  the  tube  and  of  a  magnitude  which 
varies  linearly  with  distance  from  the  center  of  the  arc.  The  assumptions 
of  a  linear  variation  of  vorticity  and  local  tangency  would  seem  to  be 
admissible  approximations  and  consistent  with  the  prior  assumptions 

such  as  the  neglect  of  shed  vorticity  effects.  The  self-induced  velocity  Vs. 

0 

at  point  PL-  --  i.  e,  ,  the  contributions  of  the  elements  adjacent  to  A  --  may 
now  be  written  in  integral  form  by  summing  the  effects  of  circular-arc 
vortex  filaments  of  differential  corss  section.  If  it  is  also  assumed  that  the 
core  radii  are  much  less  than  the  radius  of  curvature  the  integrand  may 
be  simplified  and  the  integral  evaluated.  It  is  then  found  that 


Ks. 


1 
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r, 


tan 


+  r. 


M^r-tan  ~ 

“V  T- 


(5) 


where  f>L_f  and  are  defined  in  Figure  4.  The  self-induced  velocity  is 
directed  normal  to  the  plane  of  the  arc  of  the  tube.  The  approximate  core 
radius  of  a  given  element  may  be  assigned  on  a  rotational  basis  using  energy 
considerations,  as  discussed  later. 
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CORE 

RADIUS 


Figure  4  GEOMETRIC  RELATIONSHIPS  DEFINING  THE  SELF-INDUCED 
VELOCITY  AT  WAKE  POINT  Pt 


Determination  of  System  Parameters 

The  inputs  to  the  calculation  were  nondimensionalized  for  purposes 
of  coding,  with  lengths  made  dimensionless  by  rotor  radius,  R ,  and  velocities 
by  rotor  tip  speed,  fiR .  The  flight  conditions  of  the  aircraft  being  repre¬ 
sented  relate  to  the  computer  program  through  the  following  dimensionless 

Vf 

parameters:  advance  ratio  ,  nondimensionalized  blade  vortex  cir- 

j- »  A 

culation  ^  ,  tip-path-plane  angle,  CCT,  and  the  ratio  of  wake-vortex  core 
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radius  to  rotor  radius,  .  The  quantities  which  determine  p  are,  of  course, 

R 

generally  available.  Estimation  of  the  other  parameters  is  accomplished  in 
the  manner  discussed  below. 

Radial  Maximum  of  Blade  Vortex  Strength, 

The  circulation  about  a  translating  and  rotating  rotor  blade  at  a 
given  azimuth  position  p  is  actually  varying  in  the  radial  direction  (more 
or  less  elliptically).  It  was  felt,  though,  that  a  sufficient  approximation 
to  blade -induced  effects  would  be  obtained  by  making  the  circulation  constant 
spanwise  and  of  a  magnitude  equal  to  the  maximum  attained  on  the  actual 
blade.  It  might  be  supposed  that  blade-induced  effects  would  then  be  some¬ 
what  exaggerated  but  comparisons  of  computed  and  measured  flows  indicate 
that  the  error  introduced  in  the  root-mean- square  velocity  estimate  is  not 
large  enough  to  be  discernable  (See  Reference  3).  The  error  introduced  in 
the  local  velocity  very  near  the  blade  would,  however,  be  expected  to  be 
relatively  large. 

If  the  circulation  is  assumed  to  vary  elliptically  in  the  spanwise 
direction  and  the  advance  ratio  is  in  the  normal  helicopter  range,  the  maxi¬ 
mum  circulation,  ,  about  a  blade  is  approximated  by  (neglecting  root 
effects  and  making  small-angle  approximations  where  appropriate), 


8La(p)(l-2pstnp) 

irpflRz 


(6) 


where  LB  is  the  total  lift  acting  on  the  blade  and  p  is  the  air  density.  Then, 
assuming  that  the  total  lift  acting  on  the  rotor  is  constant  and  equal  to  the 
gross  weight  W  of  the  aircraft,  it  follows  that 


=  2jrMl-2pscnP) 


(7) 
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where  A  is  a  blade  loading  parameter,  defined  by 


tfW_ 

*  "  7T2  a/. 0  p n2 r* 

with  denoting  the  number  of  blades  on  the  rotor. 

Inclination  of  the  Tip  Path  Plane,  OCt 

The  tip -path -plane  angle  is  a  complicated  function  of  flight  con¬ 
ditions,  rotor  dynamics  and  blade  pitch  settings.  Various  methods  are 
available  for  computing  this  angle  (see,  for  example,  Reference  7),  but 
their  use  requires  iterative  or  graphical  procedures  and  the  values  of 
parameters  which  are  generally  not  well  defined. 

Reasonable  values  for  (XT  have  been  obtained  without  resorting 
to  elaborate  methods  by  assuming  the  rotor  resultant  force  to  be  normal 
to  the  tip -path -plane  and  equating  the  horizontal  component  to  the  fuselage 
drag. 

cxr  sc  —DI.  iff  (8) 

TTAtg  A 

where  is  the  fuselage  drag  coefficient,  defined  as  the  ratio  of  fuselage 

drag  to  .  Equation  (8)  was  obtained  by  neglecting  rotor  in-plane 

drag  and  fuselage  lift.  The  angle  ccr,  if  assumed  to  be  small,  must  then  be 
simply  the  ratio  of  fuselage  drag  to  rotor  thrust  (aircraft  weight).  It  was 
noted  in  Reference  3  that  values  for  o;robtained  from  Equation  (8)  were 
found  to  be  within  20  to  30  percent  of  measured  values  for  a  UH-1A  heli¬ 
copter  (Reference  8)  using  a  value  for  of  0.014  (Reference  9). 

Wake-Vortex  Core  Radius,  a 

The  nondimensionalized  core  radius  of  each  vortex  element  must 
be  specified  as  it  is  generated  by  a  blade  vortex.  This  core  size  should 
correspond  to  that  of  the  tip  vortex  formed  in  the  physical  flow  when  the 
wake  rolls  up. 
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In  order  to  obtain  an  estimate  of  the  core  size,  three  major 

simplifying  assumptions  have  been  made.  First,  all  vorticity  in  the  flow 

except  that  contained  in  the  wake  under  analysis  has  been  neglected.  The 

justification  uses  the  fact  that,  some  distance  from  the  wake,  the  flow 

before  and  after  roll-up  is  effectively  the  same.  Second,  it  is  assumed 

that  the  wake  may  be  considered  as  though  it  were  generated  by  a  wing 

of  the  same  geometry  as  tne  blade  in  steady  translational  flight.  That  is, 

'I' 

the  tip  and  root  vortices  of  a  rotor  blade  are  assumed  to  have  the  same 
core  size  as  the  tip  vortices  generated  by  a  wing  of  the  same  aspect  ratio 
in  steady  translational  flight,  and  having  the  same  loading  as  the  blade  had 
when  it  generated  the  segment  of  wake  being  analyzed.  Third,  it  is  assumed 
that  the  roll-up  of  the  wake  of  a  wing  in  steady  translational  flight  can  be 
treated  as  an  unsteady  two-dimensional  problem  rather  than  a  steady  three- 
dimensional  one.  The  only  justification  (other  than  intuitive  argument)  for 
the  latter  two  assumptions  is  that  they  do  yield  results  which  appear  to  be 
realistic. 


Consider,  then,  a  wing  of  large  aspect  ratio  in  steady  translational 
flight,  as  shown  in  Figure  5.  If  r  denotes  the  circulation  about  the  wing,  then 
at  some  position  in  the  wake  upstream  of  the  point  where  roll-up  begins 
(View  A-A,  Figure  5),  the  circulation  per  unit  length  /(c/)is  given  by 

/({/) =  zf  =  uJh  (9) 

where  tti(y)is  the  magnitude  of  the  vorticity  in  the  wake  and  h(y)  is  the  thick¬ 
ness  of  the  wake  sheet.  Given  the  loading  on  the  wing  --  i.  e.  ,  on  the  rotor 
blade  at  some  azimuth  --  f(y)  may  then  be  specified  through  Equation  (9). 


The  net  effect  of  the  blade  root  vortices  is  neglected  in  the  rotor  calcu¬ 
lations  for  reasons  given  previously  on  page  4. 
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Figure 


Note  that  /=  0  at  some  point,  say  y0  ,  where  P  is  a  maximum.  Also,  since 
r(-b)  -  r(b)  =  0,  -^corresponding  to  blade  root  and  b  to  blade  tip,  it  follows 
that  y  b 

~J  r(y)dy  *J  /( y)dy  = 

Vo 

where  denotes  the  maximum  value  of  "(y).  Replacing  variation  in  #  in 
the  wake  by  variation  in  time  in  a  plane  parallel  to  the  y-}-  plane  and  fixed 
in  the  wake";  the  sheet  ox  vorticity  would  be  observed  to  roll  up  into  a  pair 
of  vortices  whose  cores  consist  of  the  rotational  fluid  which  was  cor  ained 
in  the  sheet. 


This  process  will  be  assumed,  in  what  follows,  to  be  two  dimen¬ 
sional.  That  is,  the  flow  in  the  wake -fixed  plane  is  taken  to  be  the  same 
as  though  the  wake  were  extended  to  positive  and  negative  infinity  in  the 
/^-direction,  and  that  in  each  plane  parallel  to  that  plane  the  flows  are 
identical;  or  taken  still  another  way,  it  is  assumed  that  derivatives  with 
respect  to  x,  are  small  in  comparison  to  changes  in  y  or  ^insofar  as  roll-up 
is  concerned. 


If  the  cores  of  the  vortices  formed  by  the  roll-up  are  assumed  to 
be  of  circular  section  and  of  approximately  uniform  vorticity,  and  if  it  is 
further  assumed  that  <x>is  constant  in  the  sheet  before  roll-up  (except  at  i)0, 
where  it  changes  sign),  it  follows  that  the  core  radii  of  the  two  vortices  will 
be  equal.  Further,  because  the  moment  of  impulse  of  the  flow  is  conserved, 
the  vortices  locate  on  the  y-axis  with  their  respective  position  y  =  yz  and 
=  yR  given  by 


^'pk)l  ^Y<v)dv 

i  c  b 

yY(y)<Ly 

1  *  Vo 


(10) 
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The  situation  after  roll-up  is  as  shown  in  Figure  6. 


rm 

i 

U=c/ 


3/? 


4—  ^ 

y*  *> 


Figure  6  R EPRES ENTAT i ON  OF  THE  WAKE  AFTER  ROLL-UP 


If  viscous  effects  are  neglected,  the  kinetic  energy  of  the  flow  per 
unit  depth,  T ,  must  be  conserved.  By  deriving  expressions  for  T  for  before 
and  after  roll-up  and  equating  them,  a  relationship  for  core  radius  a  may  be 
obtained. 


From  Reference  10,  T  is  given,  in  general,  by 


r= 

where 


(11) 


(y-r?)2 +(}-?) 

bz 


zi 


V; 


y 

The  integrals  are  taken  over  the  entire  area  where  cO  is  not  zero.  If  is 
substituted  for  tD  in  Equation  (11)  and  h  is  assumed  to  be  much  less  than/?, 
the  integrands  can  be  simplified  and  the  integrals  evaluated.  It  is  then 
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found  that,  in  the  limit  of  vanishing  sheet  thickness,  1  before  roll-up  is 
given  by 


r-  - 


/f. 

4-tt 


dydt!~ 


(12) 


Similarly,  it  is  found  that  after  roll-up, 


r= 


. 


2  r- 


m 


2tt 


Cb 


(13) 


so  that  the  core  radius  is  determined  from 


*  1  zj-t  J~b  r~ 


'ill 

b 


cLydrf-—  (14) 


The  result  of  Equation  (14)  may  be  put  in  a  form  more  directly 
related  to  blade  circulation  as  follows.  Let  f  be  expanded  in  a  Fourier 
series  according  to 


r(e) 

~rZ~ 


IT 


^  Cb^  Stn  hO 


(15) 


where  cos  •  Then 

b 


y(e)~- 


dr 


t>  sin  6  dd 


oo 


7T  b  Sin  9 


y'  kcb^cos  k.9 


k  --  1 


Also,  it  is  found  from  Equations  (10)  and  (15)  that 

yR_-v±  =:  —L  r  b 

n 

'  m 


rrJ  y *<*)*«- “i 
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If  y  and  -  yL  ars  then  substituted  in  Equation  (14)  and  the  integrations 
performed,  it  is  found  that 


where  is  the  radial  distance  from  the  axis  of  revolution  to  the  blade  root. 

An  indication  of  the  order  of  magnitude  of  core  radius  predicted  by 
Equation  (16)  can  be  obtained  by  considering  the  case  of  an  elliptic  distri¬ 
bution  for  r  (i.  e.  ,  minimum  wing  induced  drag,  and  hence  maximum  core 
radius).  Supposing  r0  to  be  zero,  and  substituting  coefficients  appropriate 

77“ 

to  the  elliptic  distribution  (<af  =  y,  find  a  value  for  fa  of  0.  0862. 

This  value  agrees  with  that  given  in  Reference  11. 

Some  numerical  computations  were  performed  (Reference  3),  util¬ 
izing  Equation  (16),  to  determine  the  order  of  magnitude  of  core  size  pre¬ 
dicted  by  this  method  and  to  provide  an  indication  of  the  amount  by  which  the 
core  size  varies  with  azimuth  and  forward  speed.  Data  for  the  computations 
were  generated  by  a  blade-loads  program  that  CAL  developed  for  the  U.  S. 
Army  Aviation  Materiel  Laboratories  (these  results  are  generally  in  good 
agreement  with  experimentally  obtained  loadings).  Three  flight  conditions 
were  considered:  (1)  a  UH-1A  rotor  at  105  knots  forward  speed,  (2)  a  UH-1A 
rotor  at  30  knots  forward  speed,  and  (3)  a  5-foot-diameter  experimental 
rotor  in  hovering  flight.  The  ratio  of  rotor  radius  to  vortex  core  radius  was 
calculated  for  four  azimuth  positions  of  the  blades  in  each  case,  using  a  12- 
point  harmonic  analysis  of  circulation  distribution.  The  results  of  the  com¬ 
putations  are  tabulated  below. 
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Flight  Condition 

— 

Ratio  of  Rotor  Radius  to  Core  Radius 

o 

o 

u 

n 

V O 

o 

o 

o 

o 

CO 

II 

ip  =  270° 

105  knots  -  UH-1A 

24.  0 

13.  3 

17.  5 

21. 0 

30  knots  -  UH-1A 

18.  7 

15.  0 

14.  4 

24.  1 

Hover  -  5-Foot  Rotor 

21.  9 

21.  9 

21. 9 

21.  9 

This  method  of  predicting  vortex  core  radius  appears  to  provide 
a  reasonable  oraer  of  magnitude  for  that  quantity  --  from  about  4.  0  to  7.  5 
percent  of  rotor  radius.  Most  notable  is  the  lack  of  any  substantial  variation 
in  core  radius  either  with  forward  speed  or  azimuth  position.  The  variations 
obtained  may  be  considered  small  since  a  change  of  20  or  30  percent  in  core 
radius  produces  a  much  smaller  change  in  self -induced  velocity,  with  the 
latter  varying  as  the  logarithm  of  the  former.  For  example,  with  &  changed 
from  0.  05  to  0.  065,  the  change  of  the  self-induced  velocities  would  be  approxi¬ 
mately  5%. 

It  was,  therefore,  possible  to  incorporate  a  major  simplification 
in  the  digital  computer  progiam  by  assigning  a  representative  value  of  0.  05 
to  the  ratio  of  core  radius  to  rotor  radius,  independent  of  azimuth  position 
or  flight  conditions.  It  was  verified  (Reference  3)  that  this  simplification 
does  not  cause  an  appreciable  error. 
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RESULTS  OF  SAMPLE  COMPUTATIONS 


A  few  exploratory  calculations  were  made.  All  were  done  for  a 
two-bladed  rotor.  The  table  of  conditions  covered  is  given  below. 


Case 

No. 

Advance 

Ratio 

M 

Height -to- 
Rotor  Radius 
Ratio 

H/R 

T  ip -Path - 
Plane  Tilt, 
<XT 

degrees 

Loading 

Parameter 

A 

Number 
of  Blades 
N 'b 

Tip  Vortex 
Core  Radius  - 
to -Rotor  - 
Radius  Ratio 
a> 

\ 

0 

0.  5 

0 

0.  00506 

2 

5  x  10"5 

2 

0 

1.0 

0 

0. 00478 

2 

0.  05 

3 

0.  02 

1.0 

0 

0.  00478 

2 

0.  05 

4 

0.  02 

1.0 

0.  6 

0. 00203 

2 

0.  05 

5 

0.  05 

1.0 

1. 25 

0.  00203 

2 

0,05 

6 

0.  10 

1.0 

2.  50 

0.  00203 

2 

0.  05 

Preliminary  Remarks 

The  computing  program  used  in  this  study  allowed  the  choice  of 
100  field  points  at  which  velocities  could  be  computed.  In  addition,  of 
course,  the  velocities  at  the  end  points  of  the  vortex  elements  representing 
the  rolled-up  wake  vorticity  were  always  calculated.  The  program  is  set 
up  in  such  a  manner  that  the  solution  (vortex  end  points),  after  a  certain 
running  time,  can  be  stored  on  punched  cards  so  that  additional  calculations 
can  be  initiated  without  going  through  the  whole  starting  process.  For 
example,  in  the  interest  of  economy,  it  has  been  the  practice  to  run  the 
program  for  a  time  period  that  corresponds  to  the  replacement  of  all  the 
initialized  wake  by  elements  which  have  been  free  to  convect  in  the  developed 
velocity  field.  That  is,  if  the  wake  has  been  specified  to  consist  of  m  rev¬ 
olutions  per  blade,  the  problem  is  run  for  the  time  equivalent  of  m  revolutions. 
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At  this  point,  the  computation  is  stopped,  the  selected  field-point  coordinates 
are  inserted,  and  the  problem  is  restarted  and  run  for,  say,  one  additional 
rotor  revolution.  In  this  way,  the  computations  at  the  field  points  are  made 
only  after  an  iterated  wake  position  has  been  obtained.  The  problem  can 
then  be  restarted  again  with  new  designators  for  the  field  points  so  that  the 
velocity  field  appropriate  for  a  given  problem  can  be  expanded  to  the  extent 
necessary.  One  might,  for  example,  be  interested  in  computing  the  velocity 
distribution  in  the  space  volume  to  be  occupied  by  the  wing,  fuselage,  tail, 
and  tail  rotor. 

The  particular  problem  being  attacked  also  determines  the  amount 
of  computation  required  in  the  sense  that  the  wake  shape  should  be  stabilized 
relative  to  the  space  volume  of  interest.  In  practice,  this  means  that  calcu¬ 
lations  along  the  ground  plane,  for  example,  generally  require  a  larger  wake 
array  than  the  calculation  of  the  velocity  field  very  near  the  rotor  plane. 

The  purposes  of  the  various  calculations  given  in  Table  1  are 

Comparison  of  the  distorted  helix  model  results 
with  those  obtained  from  the  vortex  ring  model 
for  hovering  flight  (Reference  1), 

Provides  a  comparison  of  calculated  results 
with  experimental  data  (time-averaged 
velocities). 

Provide  an  illustration  of  the  effects  of  a 
change  in  actuator  loading. 

Show  the  effects  of  advance  ratio  (including  the 
effect  of  the  estimated  trim  tip-path-plane  angle). 


the  following: 

Case  1 


Case  2 


Case  3 
Case  4 


Case  4] 
Case  5  r  - 
Case  6 J 
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Most  of  the  results  are  presented  in  terms  of  nondimensional 
space  coordinates  and  velocities.  In  those  cases  where  the  interpretation 
is  most  easily  made  in  dimensional  quantities,  a  tip  speed  of  690  ft/sec. 
was  used  in  converting  to  the  dimensional  velocities. 

Comparison  of  Distorted-Helix  Solution  with  Vortex-Ring  Solution, H/R  =  0.  5 
Results  obtained  with  the  distorted-helix  model  are  compared 
with  results  obtained  from  the  vortex-ring  model  (Reference  i)  for  hovering. 
Averaged  (in  the  root-mean-square  sense)  axial  velocities  halfway  between  the 
rotor  and  the  ground  plane  are  shown  in  Figure  7  and  the  time -averaged  radial 
velocities  at  the  ground  plane  are  shown  in  Figure  8.  It  can  be  seen  that  chere 
is  reasonably  good  agreement  between  the  results  of  the  two  models  in  terms 
of  the  averaged  quantities.  In  Figure  9,  the  time  variation  (in  terms  of  the 
azimuthal  position  of  the  reference  blade)  of  the  radial  velocity  along  the 
ground  at  various  radial  positions  is  shown.  Here,  differences  in  detail 
can  be  noted  although  again  it  is  evident  that  the  mean  velocities  are  approxi¬ 
mately  the  same  (note  the  broken  scales  for  *  =  0.  6  and  1. 0).  The  differences 
in  form  of  the  velocity  time  histories  are  to  be  expected.  As  the  helix  moves 
past  a  fixed  point  on  the  ground,  and  the  blade  bound  vortex  rotates  in  the 
rotor  plane,  a  relatively  smooth  velocity  variation  is  obtained.  The  ring- 
vortex  solution,  on  the  other  hand,  shows  a  discontinuous  velocity  distri¬ 
bution  at  points  in  time  corresponding  to  the  impulsive  addition  of  a  ring 
to  the  flow.  The  characteristic  dominant  frequency  for  both  the  helix  and 
the  rings  is  the  blade  passage  frequency.  Radial  ground  velocities  at  X.  - 
0.  2,  0.  6,  and  1.  0  are  in  phase  for  the  ring  solution  because  the  vorticity 
vectors  are  always  contained  in  planes  normal  to  the  shaft  axis.  The 
distorted-helix  solution,  on  the  other  hand,  contains  the  effects  of  a  swirl 
velocity  component  induced  by  the  axial  component  of  vorticity.  The  rela¬ 
tive  phasing  between  the  velocities  obtained  from  the  ring  and  helix  solutions 
has  no  significance  because  the  impulsive  ring  has  no  absolute  time  reference. 
The  ring  solution  in  Figure  9  was  phased  so  that  new  rings  were  arbitrarily 

N 'a  77- 

introduced  when  p  -  — —  in  the  helix  solution. 


26 


BB-1665-S-3 


'yX**#9,  NONDIMENSIONAL  RADIAL  DISTANCE  FROM  THE  AXIS  OF  SYMMETRY 


Figure  7  COMPARISON  OF  RING  AND  HELIX  SOLUTIONS  FOR  THE  RADIAL 
DISTRIBUTION  OF  THE  TIME-AVERAGED  AXIAL  VELOCITIES  AT 
^  =  -0.25.  CONDITIONS:  ^  =  0,  %  =  0.5. 
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Figure  8  COMPARISON  OF  RING  AND  HELIX  SOLUTIONS  FOR  THE  TIME-AVERAGED 
GROUND  VELOCITIES  AT  y  =  -0.5.  CONDITIONS:  p  =  0,  =  0.5. 
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NONDI MENS  I ONAL  RADIAL  VELOCITY  ALONG 


i 


I 


jp  ,  AZIMUTHAL  POSITION  OF  REFERENCE  BLADE,  degrees 


Figure  9  COMPARISON  OF  HELIX  AND  RING  MODEL  TIME  HISTORIES  OF 

GROUND-PLANE  RADIAL  VELOCITIES.  CONDITIONS: jjl  =  0,  2.  =  0.6. 
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At  a  fixed  point  "near"  the  actuator,  of  course,  the  velocity- 
solution  corresponding  to  Figure  9  is  very  nearly  periodic.  The  lack  of 
periodicity  with  a  fundamental  of  one-per -revolution  seen  in  Figure  9 
probably  reflects  the  finite  extent  of  the  wake  used  in  the  calculations  and 
the  tendency  of  the  vortex  elements  to  coil  at  the  wake  termination.  -  Inclu¬ 
sion  of  diffusion  and  dissipation  effects  --  which  should  have  become  signif¬ 
icant  because  the  vorticity  has  existed  for  a  time  period  corresponding  to 
roughly  60  chord  lengths  of  blade  travel  per  revolution  --  would  probably 
alleviate  this  situation. 

Comparison  of  Calculated  Hovering  Results  with  Experimental  Data 
(Time-Averaged  Velocities) 

It  was  shown  in  Reference  3  that  the  time-averaged  velocity  dis¬ 
tribution  alculated  on  the  basis  of  the  distorted- wake -helix  model  agreed 
well  with  measurements  made  in  the  wake  of  a  helicopter  rotor  in  a  tunnel 
at  an  advance  ratio  of  0.  14.  Measurements  of  wake  velocity  distributions 
for  rotors  at  low  advance  ratios  in  ground  effect  are  sparse.  In  Figure  10, 
measurements  (Reference  6)  for  a  hovering  model  rotor  one  radius  above 
the  ground  are  compared  with  computed  results.  Since  the  measurements 
were  confined  to  the  three  horizontal  planes  and  one  vertical  plane,  indi¬ 
cated  by  the  dash  lines  in  Figure  10,  the  experimental  velocity  contours 
indicated  are  not  necessarily  accurate.  The  corresponding  computed 
root-mean- square  velocity  ratios  show  the  same  general  character  as  the 
measurements  over  approximately  the  outer  half  of  the  flow  field.  The 
differences  between  the  calculated  and  measured  results  incrase  toward 
the  centerline  of  the  rotor.  It  is  speculated  that  this  difference  might  be 
due  to  the  neglect  of  the  blade  inboard  vorticity  in  these  calculations. 
Similarly,  the  measurements  in  this  region  should  be  treated  with  some 
caution  because  the  velocities  were  deduced  from  pressure  readings  and 
some  question  arises  with  respect  to  the  assignment  of  a  meaning  to  a 
"static  pressure"  reading  in  a  confused  velocity  field. 
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(a)  CALCULATED  VELOCITY  CONTOUR  MAP  FOR  A  ROTOR  1.0  RADIUS  ABOVE  GROUND 


(b)  MEASURED  VELOCITY  CONTOUR  MAP  FOR  A  ROTOR  1.0  RADIUS  ABOVE  GROUND  (REFERENCE  6) 

Figure  10  CALCULATED  AND  MEASURED  TIME-AVERAGED  VELOCITY  RATIOS.  /6  =  0  ^  =  1.0 

K 
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In  Figure  II,  the  measured  and  calculated  root-mean-square 
velocity  ratio  profiles  along  a  line  normal  to  the  ground  at  a  distance  one 
and  one-half  radii  from  the  shaft  centerline  are  compared.  Qualitative 
agreement  is  noted  except  very  close  to  the  ground  plane  where  boundary 
layer  effects  (neglected  in  the  calculation)  could  be  expected  to  be  important. 

The  method  of  measurement  employed  in  the  test  (Reference  6) 
precluded  the  possibility  of  determining  the  time  variations  of  the  velocities 
under  the  hovering  model.  The  calculated  velocities  vary  with  time  (azimuthal 
position  of  reference  blade)  and  some  idea  of  the  variations  can  be  obtained 
from  Figures  12  and  13.  Values  of  the  velocity  components  at  points  on  a 
line  parallel  to  the  x-axis  and  halfway  between  the  rotor  and  the  ground 
(Figure  12)  and  along  the  ground  (Figure  13)  are  compared.  In  particular, 
the  time -averaged  values  are  shown  compared  with  the  values  computed 
when  the  reference  blade  is  at  -  0  degrees  and  when  the  reference  blade 
is  at  ^  -  90  degrees  (see  Figure  1  for  coordinate  system).  It  is  evident 
from  the  figures  that  the  computations  indicate  sizable  time  variations. 

Effects  of  Actuator  Loading;  fj.  =  0.02,  H/R  =  1 

The  calculated  effect  of  a  change  in  actuator  loading  on  the  geometry 
oi  the  wake  vorticity  can  be  seen  in  Figure  14  where  the  traces  cf  the  tip 
trailing  \  rtex  cores  in  the  k-%  plane  are  shown  at  the  time  when  the  refer¬ 
ence  blade  vortex  (Blade  No.  1)  is  at  i?  =  0.  The  cores  for  the  more  highly 
loaded  case  (  /\  =  0.  00478)  e;  Mbit  a  wider  separation  along  a  trajectory 
than  do  those  for  che  lightly  loaded  case  (A  =  0,00203).  Further,  the  esti¬ 
mated  trajectory  is  less  inclined  relative  to  the  rotor  for  the  more  highly 
loaded  case.  These  differences  arise,  of  course,  from  the  greater  induced 
velocity  (and  the  greater  energy)  associated  with  the  vorticity  of  the  more 
highly  loaded  blades. 


*  The  small  difference  in  tip-path-plane  angles  (0.  6°)  for  the  two  cases  is 
believed  to  have  a  negligible  effect. 
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NONDIMENSI ONAL  LONGITUDINAL  DISTANCE  FROM  ROTATION  AXIS 

Figure  12  COMPARISON  OF  CALCULATED  TIME-AVERAGED  VELOCITY  COMPONENTS 
WITH  VELOCITY  COMPONENTS  CALCULATED  WHEN  THE  REFERENCE 
BLADE  IS  AT  ^  =  0  OR  <P  =  90°.  (HELIX  MODEL). 

CONDITIONS:  /i  =  Q,  =  1.0;  ^  =  -0.5,  y  =  0. 
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NONDIMEK'S!  ONAL  LATERAL  VELOCITY  Vv  ,  NONDIMENS IONAL  LONGITUDINAL  VELOCITY 


f  =  0°  x 

i  V 


TIME-AVERAGED 


X  ,  NONDIMENSIONAL  LONGITUDINAL  DISTANCE  FROM  ROTATION  AXIS 

Figure  13  COMPARISON  OF  CALCULATED  T I  ME- AVERAGED  VELOCITY  COMPONENTS 
WITH  VELOCITY  COMPONENTS  CALCULATED  WHEN  THE  REFERENCE 
BLADE  I  S  AT  ^  =  0  AND  =  90°.  (HELIX  MODEL). 

CONDITIONS:  //  =  0,  £  =  1 .0;  $  =  -i .0,  y  =  0. 
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O  BLADE  NO.  1  VORTEX  CORES,  A  =  0.00203! 

O  BLADE  NO.  2  VORTEX  CORES,  *  =  0.00203; 
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Figure  \k  CALCULATED  WAKE  CORE  LOCATIONS  IN  THE  LONGITUDINAL  PLANE  FOR 
TWO  LOADING  CONDITIONS;  /x  -  0.02,  £  =  1.0,  ^REF<  BU[)E  =  0. 
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The  reference  blade  bound  vortex  is  located  at  f  -  0°  at  the  instant 
shown  in  Figure  14.  It  can  be  seen  that  the  last  vortex  generated  by  Blade 
No.  1  is  above  the  tip  path  at  V  =  180°  --  that  is,  the  vortex  tube  is  passing 
above  Blade  No.  2  which  is  at  &  -  180°.  It  might  be  concluded  that  at  some 
time  during  the  next  half  revolution.  Blade  No.  2  would  cut  through  the  vortex 
tube  generated  by  Blade  No.  1.  This  does  not,  however,  appear  to  be  the 
case.  The  tip  vortex  from  Blade  No.  1  slides  outboard  relative  to  Blade 
No.  2  as  Blade  No.  2  rotates.  At  the  time  at  which  the  core  is  in  the  tip- 
path-plane,  it  is  also  outboard  of  the  blade  tip  and  no  intersection  occurs. 

The  miss -distance  is  small  and,  therefore,  no  definite  conclusion  can  be 
drawn  with  respect  to  the  behavior  expected  of  the  physical  system  --  the 
position  errors  introduced  by  the  idealizations  introduced  are  probably  of 
the  same  order  as  the  miss -distance. 

In  Figures  15  and  16,  the  velocity  time  histories  at  two  points  in 
the  plane  of  symmetry  and  near  the  location  that  might  be  occupied  by  a 
horizontal  stabilizer  and/or  a  tail  rotor  are  compared  for  the  two  loading 
cases.  One  of  the  field  points  (*  =  1.0,  </  =  0.  0,  ^  =  -0.  25)  is  outside  the 
slipstream  envelope  (see  Figure  14)  and  the  velocities  do  not  differ  greatly 
for  the  two  loadings  (Figure  15).  Apparently  the  differences  in  spacings  and 
vortex  strengths  for  the  two  cases  introduce  compensating  effects.  Another 
interesting  feature  of  Figure  15  is  that  both  the  %  and  ^  velocity  components 
experience  reversals  at  blade-passage  frequency.  One  could  anticipate,  ther 
fore,  that  a  tail  plane  or  tail  rotor  would  experience  shaking  forces  and  these 
would  be  transmitted  to  the  fuselage. 

Velocity  component  time  histories  at  a  field  point  inside  the  slip¬ 
stream  are  shown  in  Figure  16.  This  location  (*  =  0.  8,  =  0.  0,  %  -  -0.  25) 

might  also  be  a  candidate  location  for  a  stabilizing  surface  and/or  a  tail 
rotor.  The  calculations  indicate  that  the  relative  amplitudes  of  the  velocity 
fluctuations  at  this  location  are  nearly  proportional  to  the  loading  parameter, 
A  .  The  velocity  component  normal  to  the  reference  rotor  plane,  V^,  is,  of 
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course,  higher  for  the  higher  loading  parameter  and  this  velocity  compo¬ 
nent  dominates  in  the  slipstream.  Again,  sizable  velocity  fluctuations 
with  respect  to  time  are  indicated. 

Another  characteristic  of  the  slipstream  velocity  is  shown  in 
Figure  17  for  the  highly  loaded  case  (A  =  0.  00478).  Figure  17  shows  the 
difference  in  the  resultant  velocities  (AD.R-\fvx2-l-V^z  )  that  might  be  experi¬ 
enced  at  the  two  tips  of  a  stabilizing  surface  ( %  =  0.  8,  y  =+0.15,  ^  =  -0.  25). 
The  velocity  fluctuations  correspond  to  differential  aerodynamic  angles  of 
attack  of  approximately  7°.  Therefore,  the  oscillating  forces  would  excite 
fuselage  bending  and  the  oscillatory  couple  would  tend  to  shake  the  fuselage 
in  torsion. 

Effects  of  Advance  Ratio 

Effects  of  advance  ratio  on  the  slipstream  shape  and  the  velocity 
components  were  calculated  for  a  fixed  value  of  0.  00203  for  the  loading 
parameter,  A  .  The  estimated  tip-path-plane  trim  angle,  aT,  was  changed 
as  a  function  of  advance  ratio,  /x  .  The  pairs  (ju,  OcT  )  are  (0.  02,  0.  6°)  in 
Figure  18,  (0.  05,  1. 25°)  in  Figure  19,  and  (0.  10,  2.  5°)  in  Figure  20.  These 
plots  will  not  be  discussed  in  any  detail  --  they  are  only  intended  to  indicate 
major  features  of  the  flow  fields. 

Slipstream  Location.  Comparison  of  Figures  19a  and  20a  -~  the 
traces  in  the  plane  of  symmetry  of  the  vortex  cores  when  Blade  No.  1  is  at 
=  0°--  displays  a  slipstream  "tuck-under"  between  //  =  0.05  and  0.  10. 

Below  some  critical  advance  ratio  in  this  range  the  ground  effect  forces 
the  wake  vorticity  to  translate  upstream  near  the  wall  (Figures  18a  and 
19a);  above  this  advance  ratio,  the  free-stream  velocity  overpowers  the 
ground  effect  and  the  vorticity  is  swept  downstream.  Extrapolation  of 
wind-tunnel  data  to  free  flight  conditions  (out-of-ground  effect)  could  be 
compromised  severely  by  differences  introduced  by  this  effect  --in  particular, 
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at  low  advance  ratios,  the  forces  and  moments  measured  in  a  wind  tunnel 
could  differ  significantly  from  those  that  would  be  experienced  in  the  absence 
of  a  tunnel  floor. 

Velocities  Just  Below  the  Tip  Path  Plane.  Figures  18b,  19b  and 
19c,  20b  and  20c  show  the  calculated  time  histories  of  the  nondimensional 
velocities  immediately  below  the  tip  path  plane  =  -0.  01)  along  the  longi¬ 
tudinal  axis  ( y  =  0.  0)  for  advance  ratios  of  0.  02,  0.  05,  and  0.  10  respectively. 

* 

The  sharp  peaks  in  the  V  component  for  |^|<|  indicate  the  passage  of  the 
blade  bound  vortex;  and  the  peaks  are  of  nearly  equal  magnitude  because  of 
the  assumption  that  the  blade  bound  vortex  strength  does  not  vary  in  the 
radial  direction.  Therefore,  the  V '  values  indicated  should  be  reduced  by 
the  ratio  of  the  local  r  to  the  constant  (maximum)  value  used  in  the  calcu¬ 
lations  to  obtain  a  more  representative  Vu  radial  variation.  The  ^-compo¬ 
nents  of  the  velocity,  v  ,  show  sizable  time  variations  in  the  vicinity  of  the 
slipstream  boundary  --  a  reversal  in  sign  occurring  at  blade  passage  fre¬ 
quency  near  \x\  =  /  at  /x  =  0.  02  and  0.  05.  Interior  to  the  slipstream 
[\x\<  0.  8,  say),  Vx  changes  with  time  at  fixed  points  are  indicated  to  be 
relatively  small.  The  components  show  the  characteristic  rise  and 
fall,  at  blade  passage  frequency,  expected  near  the  rotor  plane  in  the  in¬ 
terior  of  the  slipstream.  Near  the  boundary  of  the  slipstream,  this  axial 
component  is  strongly  influenced  by  the  tip  vortex  as  well  as  the  bound 
vortex  and  large  time  variations  can  occur.  The  sign  of  the  mean  value 
changes  as  the  slipstream  boundary  is  crossed. 

Figures  1  9d  and  19e,  20d  and  20e  show  the  velocity  ti.  .e  histories 
just  below  the  rotor  (  jl  =  -0.  10)  along  the  lateral  axis  (x  =  0.  0).  General 
characteristics  are  similar  to  those  indicated  above  with  the  /,  and  4 

%  7 

velocity  components  interchanged  with  respect  to  the  effect  of  the  blade 
bound  vortex  effect. 


The  peaks  should  be  somewhat  rounded  --  the  sharpness  indicated  is  the 
result  of  using  straight  lines  in  the  plotting  routine  to  join  calculated  points, 
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Velocities  Near  Fuselage.  Velocity  time  histories  along  a  line 
representative  of  a  fuselage  location  are  shown  in  Figures  1 9f  and  19g 
and  20f  and  20g  for  /£  =  0.  05  and  0.  10.  The  lateral  velocity  component,  Vy , 
is  seen  to  be  relatively  small.  The  longitudinal,  Vx  ,  and  axial,  1/  ,  non- 
dimensional  components  are  large  and  oscillatory.  It  would  be  expected 
that  these  components  would  give  rise  to  oscillatory  fuselage  forces. 

Figures  19g  and  20g  show  the  change  in  environment  experi¬ 
enced  by  a  tail  rotor  (0.  6<’*<1.  3,  y  =  0.  0,  =  -0.  25).  It  is  apparent  that 

the  velocity  field  in  the  vicinity  of  a  tail  rotor  is  nonuniform  and,  conse¬ 
quently,  oscillatory  tail  rotor  loads  would  exist. 

Flow  Near  a  Stabilizing  Surface.  Figures  1  9h,  and  20h  illustrate 
the  velocity  variations  in  the  regions  of  the  tips  of  possible  stabilizer  posi¬ 
tions  (0.  6  <■*<().  9,  y  -  +0.  15,  £  =  -0.  25).  Figure  20h  (/<=  0.  10)  indicates 
time  variations  of  V%  and  14  at  these  locations  that  are  significant  while 
Figure  I9h  {y.  =  0.05)  exhibits  a  much  smaller  ncnuniformity.  Sidewash 
velocities,  J/  ,  are  small. 

Velocities  Along  the  Ground  Plane.  Figures  1  9i  and  19j,  and  20i 
and  20 j  show  the  velocity  components  along  the  ground  at  the  plane  of  sym¬ 
metry  {y  =  0).  The  oscillatory  component  of  K,  and  the  change  of  V  with 
jc  decrease  with  increasing  advance  ratio  in  the  region  inspected.  The 
component  is  small  along  the  line  on  which  the  calculations  were  made  and 
V  =  0  by  imposition  of  the  boundary  condition. 

o 
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CONCLUDING  REMARKS 


The  model  developed  provides  estimates  of  the  velocity  field  in 
the  vicinity  of  a  rotor  (propeller)  operating  in  ground  effect.  In  particular, 
the  time  variations  predicted  appear  to  be  reasonable  --at  least  ovei  the 
outer  half  of  the  slipstream  and  at  field  points  outside  the  slipstream.  Since 
there  are  no  known  suitable  surveys  that  show  the  time  variation  of  velocities 
in  the  vicinity  of  a  rotor,  it  was  not  possible  to  verify  this  most  important 
aspect  of  the  present  eifort  nor  to  test  the  validity  of  the  many  assumptions. 
Briefly,  these  assumptions  were  the  following: 

1.  Incompressible  fluid. 

2.  Inviscid  fluid. 

3.  Negligible  effect  from  root  vortex. 

4.  Negligible  effect  from  shed  vorticity  at  small  advance 
ratios  (/z<  0.  10), 

5.  Self -induced  velocity  determined  by  local  vortex 
curvature. 

6.  Blade  bound  vortex  constant  along  span. 

7.  Constant  initial  vortex  core  diameter. 

8.  Vortex  core  diameter  estimated  from  wing  theory. 

It  is  believed  that  the  theory  developed  is  applicable  to  a  variety 
of  problems  that  arise  from  the  interaction  of  propeller  (rotor)  flow  fields 
with  other  aerodynamic  surfaces.  It  should  be  recognized,  however,  that 
only  the  zeroth-order  effect  --  the  velocity  field  in  the  absence  of  aero¬ 
dynamic  surfaces  --  has  been  obtained.  Questions  related  to  the  inter¬ 
action  of,  say,  a  wing  surface  or  a  tail  rotor  blade  with  a  free  vortex 
have  not  been  investigated. 
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The  calculated  .'fsuits  indicate  that, 
ground  effect  forces  the  wake  to  translate1  n ^ s t ^  cam.  Consequently,  wind- 
tunnel  measurements  of  forces  and  moments  acting  on  rotors  (and  fuselages, 
stabilizing  surfaces  and  tail  rotors)  may  not  be  representative  of  out-of¬ 
ground-effect  behavior. 
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RECOMMENDATIONS 


1.  Direct  measurement  of  time -dependent  flow  quantities  in  and  near 
the  wake  of  a  propeller  (rotor)  in  ground  effect  is  recommended. 

In  particular,  velocities  in  the  flow  field  and  wind-tunnel  wall 

(or  ground)  pressure  distributions  would  be  desirable.  Comparison 

\ 

of  these  measurements  with  the  theoretical  results  would  furnish 
a  measure  of  the  validity  of  the  theory. 

2.  Vorticity  diffusion  and  dissipation  effects  should  be  incorporated 
in  the  theory. 

3.  An  effort  should  be  made  to  investigate  the  effects  of  shed  vorticity 
components  and  the  inboard  trailing  vorticity  components. 

4.  The  utility  of  the  theory  for  engineering  purposes  should  be 
established.  In  particular,  measurements  of  the  time-dependent 
pressure  distribution  on  a  wing  or  fuselage  in  a  rotor  wake 
(compound  configuration)  could  be  compared  with  estimates 
based  on  the  calculated  rotor  flow  field. 


i* 
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BLADE  NO.  2  CORE  LOCATIONS 
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Figure  18a  POSITIONS  OF  WAKE  VORTI CITY  (CORES)  IN  LONGITUDINAL 
PLANE  OF  SYMMETRY  WHEN  SLADE  NUMBER  1  IS  AT  V  =  0°. 
CONDITIONS:  ju  =  0.02,  A  =  0.(30203,  ar=  0.6°. 
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LINE  OF  LOCATION 
OF  FIXED  POINTS 
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Figure  18b  NOND 1  MEN S 1 ONAL  VELOCITIES  AT  FIXED  FIELD  POINTS  NEAR 
THE  TIP  PATH  PLANE  VS  REFERENCE  BLADE  POSITION , V  . 
CONDITIONS:  /x  =  0.02,  A  =  0.00203,  ar  =  0.6°; 

y  =  0.00,  $.  =  -0.010,  *  VARIABLE. 
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Figure  19a  POSITIONS  OF  WAKE  VORTiCITY  (CORES)  IN  LONGITUDINAL 
PLANE  OF  SYMMETRY  WHEN  BLADE  NUMBER  1  IS  AT  f  =  0°. 
CONDITIONS: ^  =  0,05,  *  =  0.00203,  ar  =  1.25°. 
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Figure  19b  NONDIMENSIONAL  VELOCITIES  AT  FIXED  FIELD  POINTS  NEAR 
THE  TIP  PATH  PLANE  VS  REFERENCE  BLADE  POSITION, ^  . 
CONDI  HONS:  M  =  0.05,  A  =  0.00203,  aT  =  1.25°; 

cf  --  0.00,  &  =  -0.010,  *  VARIABLE. 
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Figure  19e  NOND I  MENS  I ONAL  VELOCITIES  AT  FIXED  FIELD  POINTS  NEAR 
THE  TIP  PATH  PLANE  VS  REFERENCE  BLADE  POSITION, 
CONDITIONS:  ju  =  0.05,  =  0.00203,  orr  =  1.25°; 

>5  =  0,00,  f  =  -0.010,  y  VARIABLE. 


53 


BB-1 665-S-3 
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Figure  19g  NOND I MENS I  ON AL  VELOCITIES  AT  FIXED  FIELD  POINTS  NEAR  STABILIZER 
AND/OR  TAIL  ROTOR  LOCATIONS  VS  REFERENCE  BLADE  POSITION, 
CONDITIONS:  //  =  0.05,  A  =  0.00203,  ar  =  1.25°; 
y  -  -0.25,x  AND  y  VARIABLE. 
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Figure  19h 
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NONDIMENS  I ONAL  VELOCITIES  ALONG  GROUND 
CONDITIONS:  =  0.05,  %  =  0.00203,  ccT  =  1.25°; 

y  -  0.00,  %  -  ^GROUND*  #  variable. 
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NONDiMENSIONAL  VELOCITIES  ALONG  GROUND. 
CONDITIONS:/^  =  0.05,  ^  =  0. 00203,  ocT  -  1.25°; 

Lj  -  0.00,  2  =  ^GROUND’  *  VARIABLE. 
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Figure  20b  NOND I M  ENS I  ON  AL  VELOCITIES  AT  FIXED  FIELD  ,’01  NTS  NEAR 
THE  TIP  PATH  PLANE  VS  REFERENCE  BLADE  POSITION,  V'  . 
CONDITIONS:  =  0.10,  /\  =  0.00203,  ay  =  2.5°; 

Lf  =  0.00.  j.  =  -0.010,  VARIABLE. 
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Figure  20c  NOND i MENS  I ONAL  VELOCITIES  AT  FIXED  FIELD  POINTS  NEAR 
THE  TIP  PATH  PLANE  VS  REFERENCE  BLADE  POSITION,  <jJ  . 
CONDITIONS: ^  =  0.10,  A  =  0.00203,  ccT  =  2.5°: 

*  =  0.00.  |  =  -0.010,  y  VARIABLE. 
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UNE  OF  LOCATION 
OF  FIXED  POINTS 
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Figure  20e  NOND ! M ENS I ONAL  VELOCITIES  AT  FIXED  FIELD  POINTS  NEAR 
THE  TIP  PATH  PLANE  VS  REFERENCE  BLADE  POSITION,  !>  . 
CONDITIONS:  ju  =  0.10,  *  =  0.00203,  ocr  =  2.5°: 

*  =  0.00.  I  =  -0.10,  y  VARIABLE. 
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Figure  20f  NONDIMENS IONAL  VELOCITIES  AT  FIXED  FIELD  POINTS  NEAR 
FUSELAGE  LOCATION  VS  REFERENCE  BLADE  POSITION, &  . 
CONDITIONS: /a  =  0.10,  ^  =  0.00203,  ocT  =  2,5°; 

</  =  0.00,  ^  =  -0.25,  ^  VARIABLE. 
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LINES  OF  LOCATION 
OF  FIXED  POINTS 
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Figure  20g  HOND I  MENS i ON AL  VELOCiTIES  AT  FIXED  FIELD  POINTS  NEAR  STABILIZER 
AND/OR  TAIL  ROTOR  LOCATIONS  VS  REFERENCE  BLADE  POSITION, y . 
CONDITIONS:  fJ.  =  0.10,  *  =  0.00203,  oc7  =  2.5°; 

£  =  -0.25,  v  AND  <j  VARIABLE. 
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Figure  20h  NON D I  MENS  I ONAL  VELOCITIES  ALONG  GROUND, 

CONDITIONS:  ju  =  0.10,  A  =  0.00203,  =  2.5°; 

<d  ~  °*00>  £  =  ^GROUND'  J6  VARIABLE. 
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Figure  20 i 
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NONDIMENSIONAL  VELOCITIES  ALONG  GROUND. 
CONDITIONS: //  =  0.10,  *  =  0.00203,  aT  =  2.5°; 

</  =  0.00,  ^ GROUND,  *  VARIABLE. 
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APPENDIX 


OPERATIONAL  INFORMATION  FOR  THE  COMPUTING  PROGRAM 


This  program  is  written  in  FORTRAN  IV,  with  the  exception  of  subroutine 
CLEAR,  which  is  written  in  MAP.  This  routine  is  used  to  initialize  storages 
to  be  zero. 

INPUTS 


CARD  l 


NB: 

NRW: 

NA: 

NPNCH: 


NOPT: 


NTAPE; 

NPRINT: 


LNCT : 
NFPT: 
NXPT : 

NPINT: 


Number  of  blades,  NB 

Number  of  revolutions  of  wake  per  blade,  A/^ 
Number  of  azimuth  stations,  NA 
Punch  option.  If  zero,  no  cards  are 
punched  at  the  end  of  a  run.  If  not  zero, 
all  wake  point  coordinates  and  core  sizes 
at  the  final  azimuth  position  are  punched 
on  cards. 

If  zero,  the  initial  wake  configuration  is 
computed.  If  not  zero,  initial  wake 
configuration  is  read  in. 

If  not  zero,  v/ake  point  coordinates  and 
velocities  are  saved  on  utility  Tape  4. 

If  NPRINT  =  1,  coordinates  and  velocities 
for  each  wake  point  are  printed;  if 
NPRINT  =  2,  those  for  every  other  point 
are  printed;  if  3,  every  third;  etc. 

Number  of  lines  desired  per  page  of  output. 
(Not  used  in  the  study  reported  here) 

Number  of  points  off  the  wake  for  which 
velocities  are  to  be  calculated. 

Output  is  produced  at  intervals  of  NPINT 
steps;  i.  e.  ,  if  NPINT  =  1,  the  data  for 
each  azimuth  position  is  printed. 
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CARD  2 


PSIO: 

REV; 

XLAM: 

XMU: 

ALPHAT: 

RB: 

H: 


Initial  position  of  Blade  1,  ^  ^  ,  degrees. 
Number  of  revolutions  of  rotor  for  which 
calculations  are  to  be  performed,  NRV 
A 
M' 

ccT  (degrees) 

R/C  =  ratio  of  blade  radius  to  blade  chord. 
H/R 


CARD  3,  4 


GAMB: 

Al: 

A: 


Normalized  strengths  of  Blade  1;  NA  of  them. 
Core  sizes  at  Blade  1;  (NA  of  them). 

Initial  core  sizes;  (NRW)(NA)(NB)  of  them. 


Coordinates  of  points  off  the  wake  at  which  velocities  are  to  be  computed 
NXPT  points  in  all  (up  to  three  sets  of  coordinates  per  card): 

XIPT  YIPT  ZIPT 

*  y  ?- 

Initial  Wake  Configuration  -  Read  in  only  if  NOPT  is  not  zero, 

X:  (NRW)(NB)(NA)  of  them. 

Y:  (NRW)(NB)(NA)  of  them. 

Z:  (NRW)(NB)(NA)  of  them. 


A  listing  of  the  program  is  given  on  the  pages  which  follow. 


Estimated  as  Pg  (P)  =1-2  ^sin  f  in  the  absence  of  data. 
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U'if-TC  >7H WV>  !  I  S  F  ,  ^ ^ f 

c _ HikLLQLTl_  W^K!~  ynULLLT^  '■tlTTH  ^animn  fpffcts  -  main  program _ _ 

COMMON  y  (  w,>'  f  S  )  ,  Y{  7  {  M4P,  5  )  ,li(  ^4"  ,  f> )  ,  V  (  V«  "> ,  S  )  ♦  W(^4^,  s7 , 

1  r,A"A|  <  '4  %  SCVr.M.  (A  «m  (  1  A*)  t  A1  (  pi )  ,  A  ( M4^  ,  S  )  ,  p  T  ,  *  AO, 

?  VM.>  t  x-ICL  ,  XMSI.  t  fc-T  ,  m<>w,NA  .NW.NOPT,  NTAPt ,  NP&  !  NT  ,  M^VCH «  PS  !  'S 

T  x'*'l,  *|.  AM,  A  l  PH  AT  ,  DINT  ,  PS  IP,  XNA,OPSl  ,  M*1  ,  MW]  ,XNW,XNM,  TDN^, 

4  SAT,r4T,ri,r:>,PSI,TD[fyftTltT?,XJ,N,»S,JDS,lDS,loSl,XXX,YYY, 

_ M _ 7  11,  tS  T  t  I  r  ,  IrLa.V  r.1  ,  SIC?  ,  S  ir.Ttnr,r>tOFN,X^»il  ,  XNtt?,  XNU3,  TjM_, 

4  "yy  .  x  »y.  x'^’,  SIT ,  iTr,^ ,  Pa’- SO!  1  ,  SF  ,  SO! ,  SO  ,  mf  t !.  PS  ,  SF<;fi’Sl:G2  , 

7  ST', |  Nf  r  ,  yinT'!l^'  )  ,  Y!  PT{  im)  ,  7IPT  ,  VX(  l''*' )  ,vy{  IP'*)  , 

A  V7  (  ]  )  ,  ‘•VDT  ,\!A«  ,  r  AC  T'i  »  4“  A  ,  M  ,  WOP  AP  ,  S  AT  ~> ,  CAT?  ,  S  A?T  »  C  A  PT  t  THC  A  T 

•J  ,  T  •  *  S  A  T 

Oivr.'jc;  i  -m  r. 

_ f- (Ji 1 1  V  M  !  Mrr  (  \y\~  1  ,r )  _ _ 

l  CALI,  CL  ^  As>’  (  X  ,  Nr  A  "•  ) 

pi  =  ?.  l  a i. '•> 

PAP.  =  . ']7V 
TP  I  =  ? . 0*P I 

P  P  A-">  1  GOO  ,  *>  A  , '  >  •/  ,  ■•<  A  ,  \PCf  » «  m»1PT  ,  J T  A P 17  f  N° P  I M T  t  L NC  T  ,’mF  PT  ,  N  XPT »NP  !  NT  t 

I _ n'uvcp  _ 

inc"o'_F  ?!  <~) 

CALI  ^V'-'-IK  (  v| V  ^ 1 ' ;  ___  _ 

I  PINT  APf  .1  r.  ■')  wr «.•!*.•?  A 

»  P  A  n  If, Cl  ,  PSf  C>’*rV,  *1  i  ii,  X  MU  ,  A  L  °H  A  7  ,  r' P  ,  M 

n«l  nPWAT(w  .  •,) 

_ >£AP  I"V*M  II  )  ,  !  -1  _ 

->  p  a">"  1  1  t  (  '>  I  (  I  )  *  ’  “  ’  * A  5  ' 

MP  ]  =  MP 

MW  -  \v  i-1  i A 

NWl  =  NWt ! 

\]AI>  =  ‘  i-*'|A 

_ X  N  A  =  a  A _ _ _ 

0°  s  i  =  ?.''*■’!  / y  \Ta 
XNW  -  AN 

yfvlP  -  K- .1 

SAT  =  S  P  (  M  ■>!•<  A  T A'M 
CAT  =  C  P  (  A  [  ">•*  A  T  A  ^  ) 

. . S A  i  ?  =  <sAT_*«7_ _  _ _ _ 

"  *C  A  T  ?’  =  P  A  T_ ■  :•  ’? 

S  A  P  T  =  S  IN  (  ?  AL  !,<i  AT  *D  />[-) ) 

CA?T  =  <  'S|  ?.<'*  A|  oi.’ATAk  A.'. ) 

THSAT  =  l«'T 

THC  AT  =  ?/i  |*r,AT 

_ v!  f'  P  A  ^  -  S-:1^  M  ^J»_(  S“"  T  (  |  p  I  »  X‘l°  »XI  AM'S.  ->P  )  )  -XV()»»?  )  } _ 

X  m >.i s f.  V  (  y  o"*,  S  \T  )  t  f  !■  t" 

X’^wSS  =  (  xv,'if.»~  iA«>*S'.T)*SAT 
•/■'CAT  -  •■'■••  A  J  *<■  A  T 

•/OfATp  =  •<rC  A  T  /  )  i 
Cl  =  X A  T 

_ C?  =  (  V'1"-- SATt  y|  A  ') _ _  _ 

ra  =  "yi|*s'\l  *-S:'>!-  T  (  xfv  .“fvin/j  ,‘r) 

XMCL  =  M  / Y|  AM 
X ,J S I  =  \M.ICS\T /XI  A*' 

TMP  =  (■•'vnPSI 
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t  v«p  i  =  <:C'l«TtTMi't|rMpf?,i''|| 

FP.'V  =  (  TVI»-TVIP1**LP“.(  1  .r+TMPfTM-»l  )  )/OP<;r 

17  r  ao  i i ,  { ( \  c  t  7  j") .  Y- 1 ",  >!w )  ,T| = i , mo  H 

IC(-MXPT.F).  '  1  m 

PFAO  10-]  .  (  V  I  *>T  c  I  )  ,Y!I»T(  n  .  7!  PTC  n  t  !  =  l  ,.\>xprj  “  . 

103  PS  JP  =  PS  ir  +  -*o.  ''*,«  =  V 

CALL  l  O' KIT  . . 

Mf  T  =  0 

tpm“  =  ?  .  ns ;>  I/V.MP  "  . “  ” 

PST  =  PSf  it-PAi1 

psn  =  °si 

psif  =  psir*yA-)f^.''s 

I F  (  \»n p T  .  r ') .  - )  C,n  T-n  ■>  . 

2  Pl-A">  1  r  ^ (  (XT  I  ,  J  )  ,  l  =  1  , ms-'  1  )  ,  J  =  l , NM 1 

1^3  TnPMAK  1  2'~ . 

P  r  AH  1C  (  (  V(  I  ,  j  )  ,  1  =  1  ,Mvi )  ,  J  =  1  ,  \R)  _ 

RFAP  10  '■»,  (  (  M  T  .  J  )  ,  !  =  1  11  ,  J  =  1  ,M3) 

no  T.i  7 

3  0n  A  f  =■  1 ,mj 1 

XI  =  FL  m'.T  (  I-  1  )  *:}PST 

FXP  T  =  cv;>(  ath*yi  )' . 

4  on  F  J  =  ].,TJ 

PS  I  0.1  =  PSI-  +  PI  .lAr(  .1-1 
SPS^J  =  S I  V ( n  S I r  J  ) 

CP<nj  =  r ‘ i s ( p s  t  ^  j  j 

_ THT-^J  =  ATA’iPTS l^JV.1  /  PS rj_<- r  A  T  )  _ _  _ _ 

GIJ  =  (n  +  SAT-Xf  o's  J)  /^VPI-m"  " . 

c  I  J  =  S  )PT  (  ]  .  s ^^'■I’S''  J**2  1  +  XLAM*  (  PYO  J  -  1  .  )  /  {  W-'C  AT*  t  H  +  S  AT* 

1  CP<nJ)  )-X[  )/r>. <“*»■•) 

x  ( i ,  j )  =  p  t.j*c  u*fpsi  x i-rwT~J  )  f  <i  *x  n-'sr+ni  j*s at 
Y  (  I,  J)  =  -rI  K-SPU  X  I-  TWT'J  1 

_ 7  (  I  ,  J  )  =  -r  I  J»SA  T*r.)<-  (  y  j  - T ‘-t ■>  j  )  _ \  [  <-  y  »»..,;S  S  j-G  I  IfrfAT _ 

5  COMTIMLF  . ~  "  . . 

4  Cn'iTl\'Ur 

7  MPS  =  A’1'  'M  »S?  /nPS  I  f  1  .  i 

IF  { ■'■IPS  /T  ,MA 1  M"S  =  l 

on  o  j  =  i  ,\r 

_ JPS  =  M|n(NP  S+IVUIJM  )  )_/  'I  ■*  MAP  ,  M  A  ) _ _ 

"ic  (jps.p.'.py  .pc ’=  '-a . 

T  PS  1  =  JPS  _  _ _ 

on  «  i  =  i ,  'i  ■! 

IPS  =  T'->S  i  _ 

I  OS  l  =  I  PS- l 

_ I  F  (  I  P  s  I  .  f  J  .  -  )  I  o  c- 1  T  m  A _ 

GAMA  l  I,. J)  =  («,AMt'(  IPSff PS!  ,  »/?.n  . 

3  COMTP'I" 
o  C.ONT  I  MU* 

I  F  (  M  T  A  of  .  -  n  .n  j  *ji  •  t--.  j- 

WP  I  TP  (  A  )  m  i ,  '«RV  ,  N|  A  ,  ‘ "r  M  ,  V*ot  ,  MTAPF,NPR  I  NT,  l.MFT  ,  Mr'PT  i  MX~PT  ,  P  f  M  f,  " 

]  V1Vr  M,  P  F  f  ^  »  -V,  V(  AM.XMIJ,  ALP,jAT,?p,i-I 

VP  !  T  p  (  4  }  v  |  p  T  ,  v  |  r  t”  7  |  T  J 

r  A  l  L  HMI  -{  A  )  _ _ _ 

10  Of>  1  ?  1=1  ,M  ») 

on  11  1  =  1  , M  J 
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S  EG  C  1 1 J  )  =  SQRTUXU,J1-XU  +  1  ,.J  M  **2  +  ( Y(  I ,  J) -Y(  I  ♦ 1 , J )> **2  + C 7 (  I,  II- 
1 _ ?(  I»l,  J  )  \**2) _ 

11  CONTINUE 

12  CONTINUE _ _ 

13  DO  29  1=1, NW  * 

DO  28  J  =  ltNBl _ _  _ 

XXX  *  XU* J) 

_ YYY  =  Y( 1 1  J  ) _ 

III  =  Z<  I, J) 

U(ItJ)  =  0.0 _  _ 

V< It J)  =  0.0 
W(I»J)  =  0.0 _  _ 

130  00  135  IL=1,NB 

XXBR  =  X( 1,LL)»CA2T-7(1 tU ) *SA?T-THSAT _ 

YYBR  =  Yll.LLI 

_ 77BP  =  -(7  (  ItLL  ) *C  A  ?T4  X (  1  ,  L  L  )  *S  A2T«-  THC  AT  1 

XXBR  1  =  X-(2»U  »*CA?T-7(?,llT*SA?f-THSAT 
YYBR1  =  V(2 ,LL ) 

Z7BR1  =  -{  7"(2,l.L)*CA?T  +  X(2tLl.  I *<A?T*TMC AT j 

SGBR1  =  SQRT(  (  XX X-XXRR  )  **  7  H  Y Y Y- Y  YRR  )  **? ♦- (  Z7 7 -7  7 BP  )  <--*?  ) _ 

131  OP  134  IPR=1,NU 

SGRR?  =  SQRT(  (  XXX-XXRP  1  )**?+<  YYY-YYRRl  )  **?M  Z7  7-7  7t>t>l)**2) 

SEGSO  =  SFGI IRR . . 

HM1  =  SGBP  l^*?4-<;GnP2**2 
I  F  (  HM  j'.GT .  SFP5Q  )  GO  TO  13? 

HM?  =  ,?5*(  (SGBP1»SGBR?)**?-SPGSQ)*(SFGSO-1SOBP1-SGBR?)*»2)/SFGSQ 
IF  (HM2.GT.A(  IRR,LL  1**2160  to  13? 

_  GGG  =  -GAMAjrHRflU/SFGUeofm  _ 

GO  ‘TO  V?3  ~ 

132  GGG  =  -GAMA  (  I P  »  ,  LL  )  *  (  SGR°  1  4-SGRR?  )  /  (  5QBR  l  *SGR«  2*  (  (  5  GRB  1  + SGRR 2  )  **  ? 

1  "  - SFGSO  )  ) 

133  XNII1  =  (YYY-YYBP1  )  *  (  77R0-77PR  1  )-(  77  7-  77  PR1  1  * (  Y  YBR-YYBR  1) _ 

XNI.I?  =  (777-77  HR  l)#(  V  YRR-yxR-U  )-(XyX-XX3R1  ) *<  7  7 R«- 7 ’ RR 1  ) 

XNU3  =  (  YXX-Xxnu  1  )  *  (  YY'’P-YYafM  )  -  {  YYY-YYRP  1  )  -  (  XXRR-XXRR  l  ) 

Oil,/)  =  HU  ,  J)  +X'|!'1*GGG 
V(  I,  J)  =  V(  I  ,  J  )  +  XNII?  <-GGG 
W  {  T  ,  J  )  =  W(I , .))  +XNO**OGG 

_ SGBR1  =  S  G  R  R  ? _ _ _ _ _ 

XXBR  =  X  X r<  R  1  . . 

YYRP  =  YYRR  | 

7  7P.P  =  7  7PR1 

XXPR  1  —  X(  |P'-tf?f  1  |  )  *f.A?T-7  (  IPP  +  ’MJ  )  *SA?T-THSAT 

YYRR  1  =  Y(  [P-U?,LI  ) 

7  7  RIM  =  -  (  7  (  I  •-’*  +  ?,  I  l  )  * r  a  •>  t  +  y  (  j  *i  ^  t  “>  1 1  L  )  *  3  a  R  T  +  T H r.  A  T  ) 

1 34  CONTI NUf 

135  CONTINUE 

1 4  00  ?  5  l  =Tt  MR  1 

1ST  =  1 
f  NO  =  MS,' 

IfLO  =  1 _ _ _ 

IF  (L.N'B.J)  GP  TO  14 
(MO  =  I-? 

IF!  G  =  ? 

IF  ( INO.GT.'M  GO  TO  16 
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15  1ST  =  J+l 

IND  =  NW _ . 

I  FIG  =  1 

_ IF  ( IST.GT.NW)  GO  TO  18  _  _ 

16  SIG2  =  SQRT((XXX-X<TST,L)>**2  +  (YYY-YUST,m**2  +  (Z77-7(IST,L))**2J 

no  17  IR=IST,INP  _ _  _  _ 

SIG1  =  SIG2 

SIG2  =  SQRTf  (  XXX-X  (  IR  +  1 ,  L)  )  »»2  +  (  YYY-Y  (I  R_f  1  , l  H ♦»?»!  7  7  7-7  < I R»1 , L H 
1  **2) 

SEGSQ  =  SFGI  IR  t  L )**? _ _  _  _  _ 

HM 1  =  SIGl**2«-SiG2**2 
TF(HMl.GT.SFGSO)GO  TO  160 

HM?  =  ,?5*(  (S!Gl  +  5IG?)**?-SFGSO")*rSFGSO-(SIGl-SIG?J**2)  /SEGSQ 

TF(HM?.GT.A( IP  t  L )**? ) GO  TO  160 _ 

GGG  =  GAMA ( IR » L ) /SFG ( IP ♦ L ) 

_  GO  TO  161  _ _  _ _ _  _ _ 

160  GGG  =  GAMA'!  IR,lF*<SIGl+SIG2)7(  S'lGi*SIG?*U$IGUSIG?)**2-SFGSQ)  I 

161  XNU1  =  (YYY-Y  (IR  +  1>L))*(7.  (  IR  yl)  zllll  *1  t_L  li -C  ZZZ- ZIIP+1,L  \  )» _ 

.  1  ‘  ( Y ( TR ♦ L )-Y { 1 p* 1 » L i )  " 

XNU2  =  (7Z7-7UR  +  I  ,L)  >»(  X|  IRt  L)-X  (I  R»1«L  )  )-(  XXX-X It R+l , 1 )  )♦ _ 

1  (7(  IP  » L  )  -  7  (  I R  «- 1  *  1. 1  » 

_ XNU3  =  (XXX-X( IR+l,l  )  )  *  (  Y  ( I R  ,  L  I  -_YJJ_R  M  ,  L )) - I YYY-Y (I R » 1 , L ) ) » _ 

1  <XUR,L>-X(lR7f;V  fi 

U(  ItJI  =  IJ  (  I  t  J  )  +  XN1J1  *GGG  _  _  _  ... 

“VU.JI  =  Vlf.J)«-XNU2*GGG 

Wtl.J)  =  W(  It  JH-XN01*GGG _ 

17  CONTINUE 

_ GO  TO  (  in,  15)  t  I  FIG _ _  _ _ _ _ 

18  IF  (L.NE.J)  GO  TO  25 

_ I_R1_=  Lzi  _  .  _  _ 

IF  {  f  .FG.  1  )  I  0  l  =  \ 

XMX  =  (  Y  C  I  R  1  ,1  )  -  Y  (  I  R  1  +  1.1-))  M  ?1ISH1  ,1) -7  UR  1  +  2. 1  ))  -(  V(  IR1  +  I  .  L)~ 
l  Y(  IP1+2.L  )  )*{  7  (  IR1  .U-ZIIRl  +  l  ,L  )) 

XMY  =  I  7  (  I R  1  .  L  )  -  7  (  I  R  1  + 1 . )  ) ) *( X< IR1 +1  ,1  )-X(IRl  +  ?.t  ))-(Z(IRl+l.L>- 

i‘  7 ( i r i +? t l )T* ( Y ( i r 1 » i ) -x i i rV+i t l  ) )~ 

XM_7  =  (XJ  IR1,L)-X(  IRU-l.l)  )»(YUR1»1  .L)-Y(  IRl  +  2.m-(X(IRl  +  ).L)- 
l  X  (  I  R  1  +  '? ,  L  )  )  *  {  v  H  R  \  1 1.  j-Y  (  I  Ri  +  1 ,  L  )  ) 

S  T G4  =  SRG( IR 1  +  1 . 1  ) _ 

S I G  3  =  SFG ( !R  1,1 ) 

SIG5  =  SJJRJ  { JJM  I P  1  +  2.1  )-X II R 1 ,L ))**?+( Y I  I R 1  +  2  ♦  L  j -YJ  IR1.L)  )**2  + 
l  1  7  (  IP  1  +  2 » H  -  7  (  I  R  l',  l )  )  * *  ?') 

_ DEN  =  (  SIG'a  +  SIG4-SIG5)*(S[G3»SIG4  +  SIG5)»(SIG4  +  SIG5-S1G3)*(  SIG3+ 

1  *  ' SIG5-SIG4)’ ~  "  . .  . . 

_ IF  ( DEN.FQ.O.o  )  GO  TO  _ ■ _ 

I FIDFN.LT. 0.0) WRITE (6, 1^02) I , J , SI  0 3 , S I G4 , S l G5 

1002  _F  ORMAT(  ?X43  HD  F  NOMINA  TOR  NF  G  A  T  I  VF.  FOR  _R  .COMcijJA  TJ_  ON I  =I3.3X3HJ  = 

i"  I  3 »  3X6HS I G3  =F 1 6. R  *  ^  X  6H  S?  G  4  =F 16. R » 3X6MS IG5  =F16.8  ) 

RR  =  SIG3*S]G4»SIG5/SGR T ( ARS  (  DEN)  ) _ 

SOU  =  SORT  (  (  2 .  C*RR-SI G3 )  * (  R .  0*RR+S  IG3)  ) 

IF  (  SIG3»*7.LF.SIG4»»?»SIG5»*2)  GO  TO  1 9 _ , 

SF  =  ( 2.0*RR+SQl 1 ) 7SIG3 

GO  TO  20 _ 

19  SF  =  (2.0*RP-SQI1)/SIG3 

I F  (  SF.FQ.0.0)  SF  =  1.0F-2Q  _  _  _ _ 
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?0  SOI  =  M'K  !(('*.  I  ■'•>1'  \  +  S ICA  )  ) 

_ IF_  {  MCA*  *2  JL  M_S  HM  *  *2/  S ! ;*. c  *  *  ?  )  .;  1  T  P  2.1 _ 

sg~=  (?.?**.*<?  «-sof) /sir,*  . . 

r.n  m  ?■> 

?1  SG  =  (  )  /«;[■:/. 

IF  (  SC. p M  <c  =  l.'r-?'' 

??  IP  (  l  ,F  1  >  Cl.'  in  ■»  * 

nr  =  ,  j ) * ( \ i  >‘C(  •?. c* sf_/ \ i_m> mj ) )  +  .  ?5)  »g  v-'An  ,  ji*( alog( 3.1* 

I  SC/ A  {  r ,  J")  )  ♦  .  ■>’()  )"/  (  A  .  PC*  SQM  (  v  HX**?  +  XMY*$?4-X  Y7**?  )  ) 

r.n  rn  ?u  _ _ 

*■»  CP  =  (  G  A  v  A  (  I  ♦  J  )  *  (  A 1  1G  (  P .  <"**  /  A  (  f ,  .J )  )♦,?*>)  )  /  ( 4  •  ^*"RR*  S0P  T  (  XMX**2+ 

l  XMY*<t?  +  Xv7**?J  )  *  _ _ 

UII.JI  =  U(UHX«Vilf 

_ V(  I ,  J  I  =  V(  I  ,  J)»XMY*»«t _ _ 

k (  I ♦ J  >  =  w ( Ifj) f XM7*nr 

*5  CCNTINUP  _  __  ___ 

SIC  l  =  A  T  (  X  X  *  $  *  M  Y  Y  V-**  *  4  7  7  7  *  i  .  . 

XXAP  =  -  T  «l  S  A  r  _  _______ 

yypp  =  <•  . 

7  7'VJ  =  -T|if»M 

SC^w  T  =  MT(  (  Txv-'vy-Pk  **T+  (‘YYY-YY«i?  )**?+{  7  7Z-7  7RR  )**?  ) 

2A  on  2  7  1-1  »''R  . . . 

[PS  =  (vn->fMJS+  ('-A*(  l.-l  )  )  MA)  . 

fr(Lns.P'>.'')  |  pc  =  \'A 
XXP«1  =  x(  l,|  I*CA7T-T>'S  AT 

_ YYA'S  1  r  v  (  1  ,1.  ) _  _ _ 

7  70fH  =  -xT  1  .1  )*sT?f-fMc7t" 

Y'gij]  I  r  {  YV'o_YVV)*( /7GY  1 -7  7ap  )  f  (  7  7  7-77  °<?  )*(  W4P1 -yyao  ) 
x^u?i  =  ( 7  7-m-  7  7  7  ) *(  xyosi -xy.po )  f(  xxx-xxpp )*(  Z7R3 1 -77 ap) 

YMIJ1]  =  I  X  Y'-P-XXX)  *  (  YYP  7  1  -  YYn«  )  «•  (  YYY-YYCP  )  «(  XXC^l-XXPP  ) 

Sf-Pp,7  =  S  )RT  (  (  XYY- n»l> )  J  **?♦  (  YYY-YYAP  1  )  $**>  +  (  77  7  —  7  7RP  1  j  **?  j  ' 

_____  GCG  I  =G  A  Ys  (  L  17  S  )  *  (  sen  j  i  )  fj  S3  1  »SCao?*((^GPM  fCCRP?)**2-1.0)  ) 

uil»J)  =  u c  1 1 J ) -reel *v vii i  i 

v(r,j)  =  VII ,  J)-(.GCI*XM!I?1 

w ( i ,  j )  -  i ,  j i-rccj*x*nni 

I F  (  l  .rn.  I  ).|.  .ro.jjr.n  tp  ?a r  _ _ _ 

USIP^  =  r  |  n  a  r  { |_  p  s—  j  >  *  r  o  s  I 

_ 2  r^M  =  SIM  (>  S  I P  *  ) _  _ _ 

“rn$oci  =  r c s (  p s  i p*k  j 

PMMP  =  (  YXY-riiSPS  I  i  (  yyy-S  INPS  I  )  <=*2+  7  77**?. 

IP  vn?  +  s  j  C1**?.GT.  I  .M  ''r>  TD  ?5d  . 

p.MH  =  Cff’  T  (  P  Mh  ) 

up  =  .?P*f(SI(,l^PMH)**?-l.7')*(1.3-{SIG'l-PMH)**?)  . .  . 

_ IP  ( i  I^»r>  ;j>  ,  GT  ♦  1  .  M  GH  T  r  p^p _ 

. HH  =  C7VTIH7)’ 

XI IT  =  X  X X *  (  C  I M  M  * <•  ■>  +  S  I  VPSI**?/{HH*RP)  )  - YYY*S  I  NpS  I  ♦CO^PSf  *  (  1  ,Q/_ 

1  (HU*VP)_1 .  P) 

YHT  =  YYY-:  (  S  I  M’S  I  *<■  7  +  f  ^  I  **  ?J  ( HH*  »P  )  )  -XX  X*S  I  N»S  I  fef^SoS  I  *  II  ./V 
|  (n|i*£..)-l  .r) 

_ _ 7  HT  =  77  7  /  H'MfrP  °  ) _  _ _ _ 

XAlln"  r  -  Yl)  T  *  7*(  "j”,  I  )  *7.7tVy(T,  L  )  * 

Y  N|l  I?  =  -  7 -IT*  v  I  I  ,L  )  t  YH  T  *  7  (  1,1) 

XNin  =-v|'T*v(  1  ,  L)  tYaT^(  1  ,1  ) 

S  IG?  =  SCPT  ((  XHT-YI  1  ,L  ))**?♦{  YHT-YU  f  I  ))***♦(  ?WW  (  J  ,  |  »  )**2) 
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f.l)  TO  ?<,o 

?‘>H.  XMfll  _=  -YYp-7n,L)*7_//*V(l,i  )  _ . _ _ 

'Vnu>'"-  “-77?*x(  i',rT+Vv'y*7Ti  ,7") . . 

X Kj I J =  -XXX*Y( i ,1 ) tYYY*X( 1 ,1 ) 

SIC,?  =  (  Xv  x-x  (  1  1 1  )  )  **?*{ YYY-Y?  1  ft  1  )**?■*■(  7  77-/1  lf|.  I  1**2) 

259  GGG  =  C-.A"*1'  f  Lf’SJ  *{  *  I  '1  H  sir,?)  /(  c  fGl  *S  !';?*(  (  S!GJ  f  sin?)  **2-I  .'ll ) 
iij  I,  j)  =  'Id  , J) ♦x-m"1*':gg 

_ V(  I ,  J )  =  v( !  tj)»yvn7»r.r,r.  _ _ 

w ( 1 1  j  i  =  vm  f  ♦  j)“i-xf4'h*f;(.rr . . . 

GO  TO  27 

260  W  (  T  ,  1 )  =  !a  (  I  »  J  )  - r-  A  "•'*’  ( I  o c  )  * c  w p 

27  C.r\'TI\'UT 

U(  I,J)  =  U(  1  *  J  >  +  x  v»r  i. 

_ W  (  I ,  J  )  =  »:  (  I  »J  )-Y  ''_Sl  _ _ 

28  CnNTIM'IF 

29  CONTTMIJF  _  _ 

80  IFtNr.T.NF.DGO  TO  ’l 

IF(NXPT.*ir  .  )  C  A  L  L  7LCTY  _  _ _ _ 

CALL  OUTPUT 

81  NCT  =  -irifl  _ . _ 

I  F  (  NT  T  .  r,c  ,t  i>  ffff  )  r.r  j  ' 

IF(MTA?r.t/).  ;)  OP  TP  _ _ _ 

W  rM  T  r  (  4  ) 0  9  I  »  X  «  Y  |_  A  ^  ♦  A  L  PH  A  T  *  M  n  ♦  NR  «  »  N  A  i  N  W 

WRIT- (4)  (  (  v{  J  ,  J  )  ,v(  J  ,  j)  ,u  I  t  j)  ,ij(  I  ,  J)  ,  V(  I  ,  J)  ,W(  ItJKGAMMNJ),  _ 

l  A  (  I  ,  J  )  ,  I  =  1  ,  M  W 1  )  ♦  J  =  1  ?  N  a  1  ) 

_ WQ  I  T  r  (  A  )  V<,V  v,\/  7 _  _  _  _ _ 

'8ino<;i  =  d<T|*.)»'S!  “  ' . 

mpr  =  upr* i  _  _ _ 

IC  (MPS.GT.'iA)  -v,p S  =  1  . . 

I  r(  P's  I  .L  f  ,PR|F  1  r-0  TO  82  _ _ 

f  F  (  \'T  APR  .  NC .  ^  )  Fv:'>  c  1 1  F  4 

_ ic(^,n  n-i,co.~)  r.r  to  ^ _ _ _ _ 

PUNCH  l''  '4 

lO'AA  rnoM  at  (  74~i7  7  (•*  I  J^r-PTI  •  :•  U'i  VO0 TIC  I  TV  r  M_rML  AT  I  HNS - HARVFY 

18FL In  77  ) 

PUNCH  l^l  t  (  (  A(  I  ,  J  )  ,  1  =  1  ,NW)  ,  J=1  ,N8)  _ 

PUNCH  K~8f  (  (  Y(  I  ,  J  )  ,  1  =  1  ,*.V1  ) , J=1 ,N«) 

_ PIIVCU  ]_2£_2^_n  YU_!  J  )  ,  T  =  !  ,M  MiJ=l,Mt') _ 

PUNCH  KC8,  (  (  7(  J  ,T)  ,  1  =  1  ,n'!  )  i  J  =  1  ,N«) 

r,0  T'l  1  _  _ _ 

32  on  8«!  J  =  1,\M  "  . . 

T  X  1  =  X  (  1  f  J  )  _ _ _ _ _ 

TY]  =  Y(  |  ,  J) 

_ _ T  7  1  =  Ml,,l) _ 

0°  -»8  l  =  r,v'.’i 

TX?:TV]  _  .... _ _ _ _ 

TY2  =  TY1 

T  7  7  =  T  7  ]  _  _ _ 

TXt  =  X(  I  ,  J) 

_ T  Y  1  =  Y  (_I  J_J] _ _ 

"  T  7  1  =  7~(  I  ,  J  ) 

X(I,J)  =  TV7  +  XI  AM*U(  J_1  ,  J|ir'!)<;|  _  _  _ 

YIIJI  =  TY’  +  XI  AW*W(  1-1  ,  JJ^PSJ  . . " 

7  (  1  *  J  )  =  r7.>  +  XLA^'Vi(  i-i  ,  J)C8P<;! 
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IF(Y(  ItJ)*SATW(  f  •.!!*(  Al»|i.r.T.A(  1  f.l)  tr.n  TO  "*3 

X  (  Li  J  »  =  (  X»  I  t  J)*r.U-7(  I  ,J)*SAT)»CAT-(H-A(  i  ,J))*SAT _ 

7  (  I  ,  J1  =  -(  IM-aVi  ,.n  f*f  AT*(  X(T,J)  *c At'- 7  (  I  »  J ) * S AT) *S  AT ) 

cn^lT  1  NUT-  * 

x.l  =  H  PATU-1  )  *TPWn 

X(  1,.|)  =  <'<.<;(  i*s  MX.j)  _  __ 

Y{  1,  1)  =  S  INC’S!  ♦XJ) 

M_1  ,_JJ  _s_ZtT _ _ _ 

38  CONTlNilP" 

no  j  =  i  ,0‘U 

JPS  =  vci)(Nr>S+(MA‘(  J-l  )  )/t:n  +  NAn,NA) 

I  F  (  J  P  8  .  r  (,• .  m  )  JOS  =  NA 

J  PS  1  =  JPS-1 

in  jpsi  .i  ).n  jpg j_  =  ^a _ 

SFG1  =  S F  ' > ( 1 f  J ) 
r.AMj  =  r. a m a 1 7 ,  j  ) 

TA1  =  AM  f  J) 

36  nr)  *  7  I  =  ? , •’) w  _ 

GAM?  =  GAM  1  "  '  . . 

G  AM 1  =  G  A  M  A  (  )  |  J  ) _ _ _ 

GAMA  (  I  ,  J  )  =  GAM? . 

SPG?  =  SFG1  _  _  _ 

SFG 1  =  Sr G ( I ,  J ) 

SPG{  I  ,  J  )  =  S''pr  (  (  y  (  |  ,  J)-X(  J  +  t  ,  J)  }  **?  +  (  v  (  T  ,  |.) -VC  I  +1  t  J)  )*»2»m  ItJ)- 
1  Ml+lrJ))*- ?) 

TA?  =  T  A  1 _ _ _ _ 

TA1  =  A()  ,  J)  . "  . . 

A  (  I  *  J  )  =  ta->*CP.t|  sr  ;?/SF(.(  I  ,  J)  f  _ _ 

37  rriMT  I MUp 

S  FG  (  1  ,  J  )  =  SM-J  r  ((  X(  1  f  .11- X(?,J  ))**?+(  Y(  1  tJ)-Y(?t.l))*»2M7(  1JI- 

1  M?, J >)**?)  '  "  . . 

_ _ M  1  »  J  )  =  A  1  (  )  O  s  ) _ 

•GAMA  {  U  J  )  =  (f.AVIM  J"S)  fGAMu  (  jdsI)  ) /V.'' 

3  8  C^NTIMijr 

Gl?  TO  1  •»  .  " 

FAjn 
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MBFTC  /.7tr>T  LJRT.RFF 

c _ hfl  i r rn> t c «■*  h\<i-  vn p r  i c  i T v_ with  ground  fff rcrs  -  surroutinf  idout 

’  SlJRRntJTIiM-  ! **>n»jT 

X(  14%*)  f  Y(‘*4%  ,/<  14r.t.S|  »U(140»t>)fV(140fe;)fW(34r'»5)t 

l  fM  'M  n  4  %  4  )  »  R  F  G  H  4  %  5  )  ,  G  A  M  R  H  1  0  Cl  ♦  M  (  %  ♦  A  fl  4  0  *5  )  ♦’p’l  i  R  A  0 , 

7  V'*^,v^|  ,XMR|  ,  \;H  ,  M'-'K  *  N&  »  NW  t  NOPT  t  A  Pc ,  NPR I  NT  » NHVCH ,  PS  1 0  ♦ 

1  X%<,  * L  M f  A L  D H  AT  ,  f>  f  m t  *  a S  f  f'V X NA ,  0  P  S  [/NR  1  ♦  N  W 1  t  XNW  ,  MIX ,  TP\R , 

_ 4 _ SAT,f  U,r  1  ,%%  PS  1  ,TPI, XI tTl »T7tXJtNPStJPSt IPS,  )P$lt XXX tYYV. 

""*>  7  7  7  ,  1ST,  I  NT  1 1  rl  r.,RTGl,S!G?*SIG3fGGG,OFNf  XNU1  ,  XNIJ2 ,  XND%  T  R  1 »' 

4  X“X,XMV,  XM7,RIG4,  RIG%P%  ROM  ,  SFt<;oi,Sf;,nF,LPR,SFGl»SFG?, 

7  sum, i  Nr r ,  vi pt  ( ]<'<') ,  Yf  pt( i^o) , n pt(  n% ,  vx( T"o)  tW(  in?j , 

R  V  7  (  l''")  ,\*YPT,NA«- ,r- Af  r*  ,  P3 ,  H ,  A‘> ,  R  AT  ?  ,  f.  AT  ?  ,  c \p  t  fjt  a?t THC  AT 

■•?  ,TM^U 

_ ]_  WP  I  T  F  (  4 ,  1  ~~  "  >  M0  ,  ’■>  ♦  v  A  ,«'‘SI%,,SJ_F  .  XI  A**,  Y-“i|t  Al.pHATtHltH _ 

100?  "FOR^Arc'lHl  ,4  )V  rf”  »{>  r'l-ii  WA*-  T/TiT-ic  !lY"’>R'ir  MM  / /45X3  1  HN'JMRF  R 

1  OF  KLAOFR  =111  /45X11MMlMnru  op  »F\/ni  JTIOWR  OF  WAKF  = 

?  Ill  /45V  i  ]  ’  I  J;  IMfW-  0  or  A7!‘.'I,ITH  RT  AT  I  1NR  =111 
3  /4*X,  ?'‘PR!  (  IMI  TIM  ) 

3  P11.T,PH  OFf.PCT<:  /AtX'MHPRl  {  «=  I  M  A  I.  )  =  F1].M . 

4  «M  OFOCFFR  /45Y?3H1  =  f:  1  ;»  .  4  / 't  5  V  ?  3HMI J 

. .  =FV%b  74  4x'?V|7m  dma't  '  ‘  VF“nT3fHH  OFGR  F F.S "7" 

4  -*4Y?  amR  /  \  =r  1  1  .  V 

7'  41011M  =r i l . o )  '  ~~ 

?  npR  r  vm  .  r  /r|  ■)/.  T  (  ma  ) 
o  s  =  r>  #  a 

W p  1  T  r  (  .A  t  1  ^A1  )  ..  _  _ _ _ _ _ 

",Oof  fiV6va’t(  // i'  mhVr  I ",  15X??H  RTKFMTH  "'F  »M  AOF  {  ,  1  4  <  ? "'  H  f.  0  P  F  Rt7f  AT 
1‘M.ArP  1  ) 

3  nn  4 

WP  I  T  r  (  h  ,  J  '^lOS.f.A^lUl.Al  (I) 

10.''?  F  on  MM  (  nr",  ST'O.Sf  !C.M 
PS  =  P_R  +  S_ 

4  F  O-V T  I  mmc 

'<  FTlj'A’ 


J 


SIRFTC  770(1  TP  I  IST,?FF 

C _ HFL  I  f.  OPT  £  3  W4*r-  VOPTICITY  WITH  GROUND  EFFFCTS  -  SUflRTHIT  INF  OUTPUT 

SURP'ltJT  IMF  OUTPUT 

COMV.im  *(  TV.',  5)  ,v( -i4^,5}  ,  Z  {  340, S  I  ,U(340,S)  ,  V  (340 ,  R  I  ,  <-l  (  *4"  ,  5 )  , 

1  r.AWA(36f  ,  5)  ,SFT,  (340.5)  ,GAMR1  (10°)  ,~A  1  {  n~-)  ,  A  (  3  4%  A  )  ,  P  iVrAO,  " 

■>  VMP,  X-ir.L  , XMSL  ,  OR,  NPW,NA,NTW,NOPT,MTAPF,  '!'JRINT,MOvr.H,DS  to, 

3  WU,  X|.  AM,  At  PHAT,  PINT,  PS!  F  ,XA'4,»'PSf  ,  1  ,NWl  ,  X-T-f  ,  XN«,  TPN°, 

4  _ SAT, '“AT,  Cl  ,f 7,  PSI  ,TI>1  ,X?  ,  T1  ,  T?  ,  X  J  ,  OPS, J»S , !»$, I PSI ,  XXX  ,YYY  , 

5  7  7  7  ,  I  ST, IMD, IFLO, SIG1  ,  SI  G2  » SI  G%GGO  ,  OFO ,  XVUl  ,XN>»?,  X0'J3,  PI, 

6  Xwy  ,XMY,X«%  SIG4,  S!f,S,  P®,  SOI  1.,  SF  ,  S'?  I  ,SO,R%.l  PS  ,  S  Ff»l  ,SFG2* 

7  SO" ,  L'ir  T  ,  X I  PT  ( I0n  )  ,  Y I  PT  (  no )  ,  7 1  PT  {  no  )  ,  VX  (  1 )■;  VYf'loVj  ; 

■i  V  7  (  n-*)  ,MXPT,NAR,FAr.TR  ,  PR ,  H,_WO  R  AP  ,  S  A  T  ?  ,  f  AT?  ,  S  A2T  ,  C A2T ,  THC  AT 

9  , TMSAT 

1  II  IMF  =  r  _  _ _ . _ 


l  G  A  MA  (  T  ,  J  )  ,  A  (  I  ,  J  )  ,  T  =  1  ,  MW  1 ,  OPF  I  NT  ) 

IQOR  FiJPM  AT(  4YRHSTAT.  ,  1AX  1MX,  14X1HY,  1  4XIH7,  ]  4  X,?OVX,  13X?0VY,  nX?HV7, 

1  1  PXRHSTR FOGTU »  6X9HC0P. F  SI  ?/ "  /  (  I  «,*F  IS.  5  7F1S  .  R  T  ) 

ILIOF  =  1 L  INF  +■  NWl_/MAXr(MPPlN'T,n+R 
Ic(  1L  INF7g7.L*NCt7  I  L  INF  =  O 

s  r.ri^'T  in1  or _ 

TF(-\!XPT.r;o.n)  on  to  r  ”  . . " 

T  F  (  II  IOr  .  r  jn-I-'P  ITF(  4,  n^.-  )NR.,MA,NPW,XL  4%  <M»J,  A|  PHAT,OPSIO,PSIO 
WR  ITt  ( nr  a)  C  x  IPT(  I) ,  YIPTIU  ,ZIPT(I  )  ,VX(  I  )  ,VY(  I  }  ,V  %  \  )  ,  1  =  1  ,‘OXPT) 
10-04  FOP^ATI  /  f>  RX^OMV U"H  ITIFS  AT  HTHFR  f^TOTS  / 7 1 0  X 1 HX ,  1  4  X  1  Hv  ,  1 4  X  I H  7  , 
I  LAYp'lVX,  nxPMVY,  1?XP*<V7  /(c?6.5,RFlp».r)  )  ) 

4  Rf-T1K,|>I _ 

r  mo  . . . 
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MRFTC  7  7  VI.  C T  I  I S  T  t  p  FT 

c _ Hf[rcnpTi-y  wakf  viiRTiriTv  with  ground  f  f  f  f  r.  r  s  -  siirpoijtinf  vlc.ty 

'su^Thiiit  INF  vir'fY  *  ‘  ■" 

rOVMMN  *{  V.%  S)  ,  Yf  R4%  S)  ,  7  n40,S)  ,  l)  H4r  ,  5  |  t_v  (140 , 5_!  ,  W  H4Q  ,  5_)j, 

i  pamai  vn,s)  ,sfo(  5) , gampi  ( no)  ,‘aT ( Km  Vm 340,*;)  ^Pit^An,"' " 
?  VMP,  XMH  ,  t  va  t  njrw  ,  MA  t  MW »  N'lpT  t  NT  APE  ♦  NPR INT  t  NDVCH  »  PS  I  0 » 

3  Xmij,  yl  AM,  A!.:’HATtP!*T,t\STF,XN\,nPSI  ,  NftT,  NW1  ,_XNW  ,  XNBT’TPNB  t 

_ 4 _ *>AT  t  C  A 1  t  f  1  «r?yP  ST  ,TpI  ,  XI ,T1 , T? ,XJ  ,NPS,  JPS,  IPS,  I  PS  1  ,  XXX  ,  YYY  , 

*>  7  77  »  1ST,  I  NO »  I  cl  G.SIGl  tSIJ2»SlG^,r.Gr,,nENtXNUl,XNU2»XNU3tIRlt 

6  XVX,XMY,XV7  ,SIG4,Sir,S,«P,SQI  l.SF,  SOI  , SG » BF , L PS ♦ S FG l ♦ SEG2 , 

7  SIJM ,  L  Mf  T  ,  X I  p  T  (1  )  ,  Y I  PT  j  1™ ) ",  7 1  Pt  H  OP)  »  VX  ( lOO  )  »VY(  lOOTT 

S  V?(  l^o)  tNXPTfNAHtFACTPfR‘i#HtwnBAR.SAT2,CAT2»SA?T,CA?Tf  THCAT 

0  ,  THS A  T  .  "  “  - 

OIMFNSIf'N  GAMP  (ICO) 

EQUI  VALFMCMGAMPl  ,GAMP  )  “ 

1  OH  7  1  =  1 » f J  X  P  T 

XXX  =  X  I P  T  {  r ")  ■■  "  . .  .  '  "  . 

YYY  =  Y  I  PT (  f ) 

7  7  7  =  7  IPT(  j  >  “  '  '  . . 

_ VX(  I)  a 

VY ( I )  =  r.* 

V7  (I)  =  2*° 

no  no  135"  Tl  =  i  '  "  . . . . . .  “ 

XXPP  =  X  (  1  , 1.  L  )_*CA?  T- y  (  1,|,  l)*SA?T-THSAT 
YYRO  a  vTl  ,ILV  '  “ 

_ 7.7RR  =  -(  U  1  tlL)*CA?T4-V()  ,LL)  *S\2T+THCAT) 

X X 0 p  1  =  X(?,Ll)*CA?T-7(?,|  L)’*SA?r-THSAT 
YYRP!  =  Y  (  ?  ,  |.  L  ) 

77RP1  =  -(7(?,Ll  )*CA?T+X(?,l  L)*SA2T<-fHCAT) . “ 

SGPR  j  -  so-;  t  (  (  XXX-XX™  )  **?  +  (  YYY-Y  YBR  )  **?+_(  7  7.Z- 7  ) 

m  no  i 74  |oj=i,nw 

_ SGRR  ?  =  _SOPT  (  (  xxx-  XXRp  1  }  **?  +  (  YYY-YYRP  1  )  **2  *  (  7  7  7 -7  ?  RO  1  )  **2  ) 

SFGS'.'  =  S‘?G(  PR  '  ’ 

HM  l  =  SO.np  l  S0RR  ?**? 

IF  (H*M  ,f.T.S  =  G«;A)  r,r,  Trj  \  ??  . . 

HM?  =  .?R*(  (SGPPl  +  SG0’"f7)*«7-SFGS0)*(SFGS0-(  SGRR1  -SGRR?  )**2  ) 7SFGS0 
IF  (Hv?.p,T.A(  IPO, LI  l**?)GP  TO  |B2 

_ GGG  =  -r.  Am  A  (  I  R»  ,  L  L  )  /  SF  G  (  I  •>»  ,  L  l  ) 

GO  TP  nR 

13?  GGG  =  -GAmm  JOR  ,u  )  *  (  SGPP  1  fSORP?)  /  (  SGRm*SORP?*(  (  SGPRH-SGPR?  )**? 

1  -Sr  GS-P  )  ) 

PI  XN!J1  =  (  YVY-YYRP  1  )  <•  (  77PR-77PR  J  )-(  ?  7  7-77RRl  )  <«(  YYRR-YYRR  1) 

XNU?  =  (  77  7-77RP  1  )*(  XXRR-XXRRl  )-(  XXX-XXRR1  )4(  77RR-77RRi  )' 

_ XNUR  =  (  XXX-XXGP  TJjM  YYRP-YYRR1)-(YYY-YYPR1)*(XXRR-XXRP1) _ 

VX(I)  r.  VX  (  I  )  +XNU1  *OGG 
VY (  I  )  =  VY{ I)+XNU?*OGO 

V 7  (  I  )  =  V7(  I  )  +  XNUl^GGG  . .  "  ' 

SGRR  1  =  S  G  R f’  ? 

XXRR  =  vyP'M  *  **"  '  “ 

YYRR  -  VYJ12 1 _ 

7  7  PR  r  7Tf,0|  " . 

XXRR  1  a  X  (  I  -  ■<  <■  ? ,  L  L  )  *C  A  ?T-7  (  I'7R+7,H  )  v  SA?  T-THS  AT 
YYRR  1  =  Y(  I»R  +  ?,  U) 

7  7  Pu  1  =  -<  7  (  IRR  +  P,l  l  )*f  A?T  +  X(  IRRf  P,U  )*SA?Tf  THCAT) 
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1  34  CONTiMlIF 

n5  rfmrnif _ 

z  on  5  j  =  i  ,'j'u 

SIG2  =  SO«T(  {  XdT(  I  ) -X  (  1  ,  J)  )**?M  Y!PT(  I  )  -  Y  (  1  »  J)  \  **'»f  (  7  1PT  (  I  )- 

1  . Z(T, 

3  DO  4  _ K  =  1 fNW 

S  ! G 1  =  s  ir.?' 

siG2  =  soph  (  xipt  ( i  )-x(K4-i  tj ))»*?»(  yi  or  (  t  )-v ( k* b^LL)j^?iiITPT ( 1 1 : 

1  7(K4l, J) )**?) 

SFGSQ  =  Ser-fK,  J  )<=*? 

H«f-  s  igT**?+s  ig?**? 

I  F  (  H  M 1  .GT.STGSn  )G-J  TO  3'" 

H«?  "=  .  ?5*(  (  SIG1  +SIG?)**?-SEGSO)*(  SFOSO-  (  STM- c  ?r. ?)$*?)  /  S«T.S0 

IF(HM2.GT.A(K»  J  )**?>GP  TO  3? _ _ _ _ 

GGG  =  GAMMK,  J)/*EG(K»J) 

GO  T_n_3 1 

30  GGG  '  -  GAM A ( K»  J  )  V("s  I G fVs  1 G?  T  /  (  p  IG 1  *S  f  r>? {  (  S  f  G 1  *  S  I  r  ■> )  =*  *?-S  rGS° )  ) 

31  *NU1  =_(  Y  IPTH  l-y(  X  t  J)  )*(7(X,J)-/{**l,J))-(MOT(J>-7(Ktjn* 

1  mK"fJ)-V(K*],J») 

XNU?  =  (7IPT(I  )-7(KtJ)  )*  (■/(«  Lj  )  -X  (  K  f  1  ,  J  ))  -  (  Xp>M  ')-  <  d  ,  J  )  )  *  _  _ 

1  (7(K,J)-7 (Kf  1,  J  )  j  ‘  '  . . . 

XNU3  =  (  V  I  PT  (  I  )  -  X  l  K  f  J )  )  *  (  v  (  ><  ,  j  )  -y  (  *  H  ,  J  )  )  -  (  Y  l  <>  i  (  I  )  -Y  [  *  ,  J  ;  )  * 

1  ”  (X(K, JJ-XIK  +  1  , J)  ) 

VX(  I)  =  V  X  (  I  )  4-  V  Ml  1 1  *■'(!, r. 
vy Cl)  '"s'  VYC  T )  +  you?*rr,r. 

V7  (  I  )  s  V 7  (  I  )  4XMM?«G„G  __  _  _  _ _ 

4  CONTINUE  ‘  "  '  . . .  . ~ 

_J  CONTINUE 

SIG1  =  SOOT (XIPTCI )**?+YTPT( I )**>+7lPT( I )**?) 

XXPR  =  -THSAT 
YYRR  =  0.0  ” 

Z7RR  =  -THCAT _ _ 

SGBp  1  =  SOHl  (  XXX-YX°f«  )**?♦(  YYY-YYP'J )  **?♦'(  7  f7-7  ?">'  )**?  ) 

_00_  6  l  =1  ,MP 

L  P$  =  v^n(M9S+(NA«(  L-J  )  ) /V-3  +  NAH.NA)' 

JFdPS.rO. LPS  =  N  A 
X’XBR]  “s  /'(I  dT*CA2T-THSAT 

YY3P  I  =  Y {  1 1 1. ) _ _ 

7  Z  30 1  =  -X( 1 ,L)*SA?t-THrAT 

XNU11  »  (YY:3p.-YYY)*(  T  /  RP  1  -  7  7  3  R  »4  (  7  7  7 -7  7  n  0  )  *  (  YY3,<  l  -YY'H  ) 

XNU21  =  I  77  33-7*7  7  )*  (  XXRRl-XXOP  )  «•(  XXX-YYOR  )*(  77t'>]  -7  Mp) 

XNU31  =  l  XXMR-XXX  >  $  <  YY«P  1  -  YYKP  J  ¥  (  VYY-YYPR  )  $(  X  X3->» -X  VRP  ) 

"SGRO  ? ""=  SORT  f  (  XXX- XV  )l*  \  )  **•?  t  (  yYY-YY  ’R  1  )**?+■(  7/7-7  7  JR  !  )**?.) 

GGG  1  =GAMR(  l  PS  )  *  (  SGR  J  i  +  s»:«  l->)  /  (  SGn,?l  vSrR»->«(  (  soroi  tSGR-3?  )*<=?-  1 .0  )  ) 

VX  {  I  )  =  VX!  T  l-GGGl  *x\|'ll  1  —  •  -  --  —  -  - . -  - 

VY  (  I  )  =  VY  (  I  T-GOG 1  *  XVM<?  1 
vzlii  =  V7  {  I  )-i.GGl*X\il'd 
XNU  1  =  -YT»T(  ()  *7  (  1  ,  L  >  +  7  M»T  (|  }*Y(  l  ,!  ) 

XNU?  =  -7IPTm*X(l,l  )  +  X[PT(I  )*7(  l  ,  i  ) 

.  XNU3..^U-y  IPT  (  I  )«v(  1  ,1  )»Y[”T(  l  )»X(  I  ,1  1 

S  f  G  “7  =  S  H’  T  (  (  X  I  P  T  (  l  i  -  X  (  1  ,|  )  (  -’TC Y  i'^>  f  (  r )  -Y  M  ,  T  n  *  *  ^  t  (  7  f  “  M'  T )  - 
1  7(1  ♦!))-*■*?)  • 

pfm  =  sp;i*s((;?*((si  ;itsi'-->)**->-i.A) 

IFfpFN.F-?.  '.n>  Op  *■'  =  .“-'l 
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gi;g  =  gam'-.i  (i  [><;>*(<;  101  ♦  sir,?) /neM 

V  x (  1  )  =  vx (  !  ) ♦  XNIJ1*GGG 

v y {  n  =  vy m ♦  x'fJu?*;Gr, . 

V/M  }  =  V 7  (  I  )  +  XI\!in*GGG 

6  rl’MT  T  Mi ir 

vx<  I )  =  VX(  n  + Y«ri 

V/(I)  =  V  7 ( 1  )-XMCL 

LSDlll  LNJlr. _ ... _ 

5  r  tiTG'm . . 
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